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Title of 
Invention 



Optimal Surface Mitigated Multiple Targeting System (OSMMTS) 



As the below named inventor, I declare that: 



This declaration is directed to the invention titled: " Optimal Surface Mitigated Multiple 
Targeting System (OSMMTS)" 

I believe that I am the original and first inventor of the subject matter which is claimed and for 
which a patent is sought; 

I have reviewed and understand the contents of the above-identified application, including the 
claims, as amended by any amendment specifically referred to above; 

I acknowledge the duty to disclose to the United States Patent and Trademark Office all 
information known to me to be material to patentability as defined in 37 CFR 1 .56, including 
for continuation-in-part applications, material information which became available between the 
filing date of the prior application and the national or PCT International filing date of the 
continuation-in-part application. 

All statements made herein of my knowledge are true, all statements made herein on 
information and belief are believed to be true, and further that these statements were made 
with the knowledge that willful false statements and the like are punishable by fine or 
imprisonment, or both, under 1 8 U.S.C. 1 001 , and may jeopardize the validity of the 
application or any patent issuing thereon. 



FULL NAME OF INVENTOR: 



Inventor: Mr. Timothy Grant Hall 


Inventor 


Signature : TGH 


Citizen of : US 
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endproc: 
This code is used in chunk 12. 



17, Target Position Report Example #1 

18. <TPR Example #1 SDU Coordinates 18) = 
{xl,yl): = (2,3); (x2,y2): = (4,2); (x3,y3): = (3,1); 

This code is used in chunk 12. 



19. Target Position Report MAPLE Calculation For Example #1 

> c:=l; (xl,yl):=(2,3);(x2 J y2):=(4 J 2);(x3,y3):=(3,l); 

c:= 1 
xl . yl :— 2, 3 

x2, y2 := 4, 2 
y ,? := 3, 1 

> osmmts(tl J t2 J t3) [1] ; 

20. A- Value Of Target Position Report For Example #1 

1 l"\%4r 3 tl 2 tl 

-ktl - ft?) c((2 tl-t2- ft?) c + ((-(50 - * - 26c 4 tl 3 fcg - 30 tl c 4 ft? 3 

6 /cl 

+ 16 c 4 tl 2 tj? 2 - 30 tl c 2 ft? - 12 c 4 t£ 3 tl + 26 c 2 tl ft? + 70 c 4 tl 2 t3 2 
+ 10 ft? 2 c 4 t2 2 + 10 ft? c 2 t2 - 8 ft? c 4 t2 3 - %4 ^ + 16 ^ 

_ 4%4? tlV2 2 + 7%4^ _ 62 c4 ^ 3 ^ _ C 6%4c 3 tlftn/ + 1Q ^ ^ tf 

7c 1 7c 1 7c 1 

4%4t3c 3 t2 2 . 4 ctA 5%4cft? %4ctl 5 c /c4c 3 ft? 3 4%4ct2 

+ ^ h 0 C 4 ft? 4 -- h — ^ ^ h ■ 



%1 ' %1 %1 %1 %1 

+ 2 % f W - 18 c 2 tj? 2 + 2 c 2 tl 2 + 46 c 4 tl 2 t3t2 + 4 c 4 ft/ ft? 2 ft? 

7c 1 

- 50 c 4 tl ft? ft? 2 + 10 c 2 ft? 2 + 22 c 4 ft/ 4 + 5 c 4 ft? 4 ) /(%1) - (50 - 26 c 4 ft? 3 V2 

- 30 tl c 4 ft? 3 + 16 c 4 tl 2 t2 2 - 30 tl c 2 ft? - 12 c 4 ft? 3 ft/ + 26 c 2 tl ft? + 70 c 4 tl 2 ft? 2 

+ 10 ft? 2 c 4 ft? 2 + 10 ft? c 2 ft? - 8 ft? c 4 ft? 3 - 62 c 4 ft/ 3 ft? + 10 ft? 3 c 4 ft? 
4°/c3c 3 tlft? 2 77c3c 3 tl 2 ft? 13 7c3c 3 tl 2 ft? 6%3c 3 tlt2t3 



I 



c /cl c /cl °/cl %1 

16 %3 ft? c 3 ft? 2 4 0 /c3ft?c 3 ft? 2 %3 ft? 2 c 3 tg . 4,^4 %3cft/ 

+ %i ~ + %r~ %i +oc L + %i 

0 °/c3 c 3 ft? 3 4 %3 c t2 2 °/c3 c 3 ft/ 3 _ 5 °/c3 c ft? , 2 , 2 
°/cl + %1 %1 °/cl * " + 

+ 46c 4 tl 2 ft? ftS + 4c 4 tl t'2 2 ft? - 50c 4 tl t2 ft? 2 + 10c 2 ft? 2 + 22c 4 ft/ 4 + 5c 4 ft? 4 )/( 

%1)) + c 2 ft/ 2 - 2 c 2 tl ftg + c 2 t^ 2 )^ 2 ') + ^ + I(« - t2) c((-(50 

- 13%4 f/ i2 ^ - 26 c 4 ft/ 3 ^ - 30 ft/ c 4 3 + 16 c 4 tl 2 ti? 2 - 30 tl c 2 ft? 

7c 1 

- 12 c 4 ft? 3 tl + 26 c 2 tl fcg + 70 c 4 ft/ 2 ft.? 2 + 10 13 2 c 4 f2 2 + 10 ft? c 2 ft? - 8 ft? c 4 ft? 3 
_ %4 13 2 c 3 tg 16 %4 ft/ c 3 ft? 2 _ 4 7c4 c 3 tl t? 2 7 °/c4 c 3 tl 2 ft? _ 4 , 

°/cl c /cl ^ 7cl %1 

- 6%4C 1"^^ + 10 tfV * + 4%4tgc 3 t. 2 +5c 4 tg4 _ 5%4c^ 

7c 1 7cl 7cl 
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By (13), we have 

2 (,/,-/„) MUhM) 

f s *2_ s *2\_/ s 2 _ s 2 A 

with a similar statement for v 21 2 , L — — ~ A'h (tl,t 2 ,t$). Let 



65 - bt ( 65 



(X* 0 ,y* 0 ) = ( ^ ,771! ( ^ - 1 ) + 61 

\mi-m 2 \nii-m 2 J 
with similar results for other forms of (#oi2/o) corresponding to (xo, yo), where 

1 2{y 2 - yi ) 

, * _ S 22 - *12 + (*3 - *l) + (»3 - 2/l) 

2 2(y !} -y 1 ) 

Note how the value of mi and 7n 2 are unaffected by consideration of the i|. 
Then 

1 2(y 2 - yi ) 2(y 2 - yi ) 

s 2 22 -sl 2 + {4-xl) + (yl-yl) {s* 2 l-s*£)-(s 2 22 -s 2 12 ) 



or 



2 " 2 (ya - Vl ) 2 (i/ 3 - Vl ) 

h * _ h , ( S 21 _ _ ( s 21 _ 5 ll) 



65 = &2 + 



( S 22 ( S 22 S 12) 

2(2/3-2/1) 



By Definition 10, this means 
Therefore, 

and by Definition 11, we have 
Likewise, 



bl-W~Mxitl.tl.tt) 

b;-h~ M^ti.ti.tt) 



xl~x 0 — 





-»T 


6 2 -6i 


mi 

- (4- 




777.1 — 7712 
-(62-61) 


_ (bi- 


mi- 


-7712 

-(6f-6i) 




777.1- 


-7712 



(mi - 7?i 2 ) (x* - x 0 ) ~ M 2 {tl,t*, t*) 



1" 



and by Definition 11, we have 



-2/0 - "i + (mi (*z*)+*0 

m . ( b i-^ ^2-61 ^ 

\ 7711 — 7712 7711 — 7772 / 

^((^-M-(^-M) 



l"!?)(vS-tt)) ~M 2 (^,«) 
r?ii 



30. The Error Likelihood Ellipse 

Definition 13 TTie Error Likelihood Ellipse fiJiii^ for the OSMMTS position (#o,2/o) w 

/ 1 / 2/0-2/0 V _ R 

with (mi - m 2 ) (xq - x 0 ) ~ M 2 (t*,t 2 ,ts) and (l — ™^ (y$ — i/ 0 ) ~ M 2 (tl.t^.tl). and where the constant 
R depends on the two M 2 {tl,t 2 ,tt) distributions. 
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since the angle of incidence a would be consistent "inbound" to, and "outbound" from, the obstruction. 
Define the Window for Reflections Between I and T to be 

^min ^ , ^ <^max -\ 

<t < (lo) 

C C 

Note that at point (x s ,y s ) we have 



a — arccos 



(x 0 - x s ) * (x r - x s ) + (y 0 - y s ) * {y r ~ Vs) 

, \J ( x o - x s) 2 + {yo - y s ) 2 \J (x r - x s ) 2 + {y r - y 8 ) 

so that 

cos (ir — 2a) — — cos 2a 

— 1 — 2 cos 2 a 

I _ ^((^o—x 3 )*(xT—x,)+(y 0 -y,)*(y,.-y s )) 2 
((x 0 -x 5 ) 2 +(i/o-: !/ s ) 2 ){(x,—x s ) 2 +(y r - Vs ) 2 ) 

((xo - Xl) - s (x r - Xl)) * ((1 - s) (x r - Xl)) 
+ {{yo - Vi) - s {y r - yi)) * ((1 - s) {y r - y,.)) 

{{XQ - Xl) - S {X r - Xl)) 2 \f ((1 - 5) {X r - Xl)) 



= 1-2 



= 7 00 

and 



+ {{yo - vi) - s {y r - yi)f )\+{{l-s) (y r - yi )f 



di («) - Y {{xi -x,_.)-s (x r - xi)) 2 + {{yi. - yi) - s {y r - y,_)) 2 



and 

d-2 («) - V {{xo -x,_.)-s (x r - xi)) 2 + {{y 0 - yi) - s (y r - y,_)) 2 
which means d 2 (s) may be viewed function of s. 



48. Policies 

The Mitigation for Reflections Method calculates 



imin _ i min fd +d 2 } and rf ™f* — \ min \di + d 2 } 

c c o<s<i c c 0<s<1 l 



subject to 

d 2 (s) = df (s) + 4 (s) - 2d! (s) d 2 (s) 7 (s) 
for each combination of SDU position, possible target position, and L and R positions for a straight line 
obstruction that would potentially corrupt the TPR through the detection of reflections (see again Figure 
10). 

Given a TPR (x 0 , y 0 ) and multiple detection times {t±, i 2 , t 3 , . . .} at SDU I, where ti is used to calculate 
the TPR, if 

^iniii ^ ^ ^ djii&x 
c — t — c 

for all i — 2,3,..., then the TPR is accepted, i.e., there is no need for a mitigation for reflections 
If, however, 

U < or < ti 

c c 

for any i — 2,3,..., then the TPR is rejected, and a new TPR is calculated using the next highest priority 
qualifying combination of SDU's. This cycle continues until a TPR is accepted or the qualifying combination 
of SDU's is depleted. 



49. MAPLE Implementation Code 

50. { Mitigations For Reflections 50 ) = 
{ Mitigation Basis Algorithm 51 ) 

{ Mitigation Extrema Algorithm 52 ) 
This code is used in chunk 12. 
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abs: 34. 

All: 16, 34, 51, 52. 

ans : 34. 

arccos: 16. 

assigned,: 34. 

betal : 16. 

beta2: 16. 

bl: 16. 

fc?: 16. 

cdf: 34. 

cni: 34. 

conjugate : 16. 

Consulting: 16, 34, 51, 52. 

convert: 51. 

Copyright: 16, 34, 51, 52. 
ctM: 16. 
(2: 51. 
<U : 51. 

d. 2: 51. 

description: 16, 34, 51, 52. 
dt: 34. 
dl: 16. 

enrf: 16, 34, 51, 52. 
eval: 16, 34. 
evalf : 34. 
el : 34. 

e, 2: 34. 
e3: 34. 
floor: 34. 
/rem: 34, 52. 
fsolve: 51. 

16. 
12: 16. 

5= 34, 51. 

global: 16, 51, 52. 

#r: 34. 

gy. 34. 
-tffc : 34. 
tra : 34. 
my: 34. 
J: 34. 

34. 
j2 : 34. 

34. 
Atf: 16. 
kOx: 16. 

/oca?: 16, 34, 51, 52. 
M: 34. 
maximize: 51. 
maxres : 52. 
maxx: 34. 
maa;?/ : 34. 



The Index 

mm: 16. 
minimize: 51. 
minres: 52. 
mtra : 34. 
m-my: 34. 
mi%: 51, 52. 
msh: 34. 
MOGEM: 34. 
ml: 16. 
m2: 16. 
nop5: 16. 
normald : 34. 
nl : 16. 
n2: 16. 
n5: 16. 
op: 16. 

options: 16, 34, 51, 52. 
osmmts: 16, 34. 
Pqi: 16, 34, 51, 52. 
proc: 16, 34, 51, 52. 
res: 52 . 

Reserved: 16, 34, 51, 52. 

K«?Ms: 16, 34, 51, 52. 

simplify: 51. 

5mm: 34. 

so/ve: 16. 

s<jrt: 16, 51. 

statevalf : 34. 

stats : 34. 

string: 51. 

swfo: 51. 

substring: 51. 

50: 51. 

si: 16. 

s2: 16. 

5,?: 16. 

then: 16, 34, 51, 52. 

thO: 16. 

io: 34, 52. 

tx: 34. 

type : 34. 

tl: 34. 

fc?: 34. 

tS: 34. 

undefined: 16, 34. 

witt : 34. 

xJ: 51. 

x.r: 51. 

x_5: 51. 

X-0: 52. 

: 51. 
xmit: 52. 
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x 2 - 2xxi + x\ + y 2 - 2yyi + y\ - s'( 

x 2 - 2xx 2 + x 2 , + y 2 - 2yy 2 + y\ — 4 ^ 

2x (x 2 - xi) + 2y (y 2 - y±) = (sf - s|) + (x 2 , - x\) + [y? 2 - y 2 ) 

2x {x 2 - Xl ) + 2y {y 2 - y L ) = ±fc^4n + (/c) 2 + (x 2 - x{) + (y 2 - y 2 ) 
Now if i/i — y 2 , then we have 

2x (x 2 - xi) = ±/c^4k + (/c) 2 + (4 - x\) 



or 



±fcJ4K + (M 2 + (x 2 l -x 2 1 ) 

X — 

2 (x 2 - x^ 

which is a vertical line. If y L ^y 2 , then we have 

2x (x 2 - x^ + 2y (y 2 - Vl ) = ±/c^Ak + (/c) 2 + (4 - x 2 ) + {y 2 2 - y() 

-/< V i'.- + U'-r + {4 - 4) + {vl - vl) 
2(2/2-2/1) 



or 



y--x 



x 2 -x L \ 
2/2 - Vl ) 



+ 



which is a line with slope — and j/-intercept ± ^ c V^ 4 *+(/ c ^+(*2^ J, i)+G'2 3/i)_ rp Q sunl marize, } et 



Li 



x — 



±Jc,yfAK+(Jc,f + {xl-xl) 



2(X2-Xl) 



y = -x( 



+ 



±/cV4 K +(/c) 2 +(xj-xj)+(yl-y^) 
2(1/2-3/1) 



2/1 = 2/2 

. 2/1 ¥■ y-i 



(7) 



is a line in x and y. Hence, as 9 varies, the set of points of intersection between Circles I and II is a line 
defined by L\ (see Figure 5). 



1H ~- 



M II (x 2 , V2) 




Figure 5: Foci of Points as $ Varies 



If a third SDU, named III, is added to the circumstance, at position (#3,2/3), with sides relabeled to 
account for the third SDU, and which is not coUinear to the positions of SDU I and II (see Figure 6), and 
subject to the axioms mentioned earlier, e.g., see Axiom 2, then another line of points of intersection L 2 , 
defined as in (7), may be expressed in terms of the same equations as before, this time for Circles I and III 
(see Figure 7). In particular, the equation for Li becomes 

±hc,yj4.K. 1 +(J 1 c.f+{xl-xl) 
2(x 2 —x\\ 



X — 



Li : 



±/icy , 4/ti+(/ic) 2 +(x|- 

2(3/2-3/1) 



-xp + tyl-yl) 



2/1 - 2/2 
2/1 ¥■ 2/2 
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%1 %1 c /cl °/cl 

+ 46c 4 tl 2 ft? ftg + 4c 4 tl t2 2 ft? - 50c 4 tl t2 13 2 + 10c 2 ft? 2 + 22c 4 tl 4 + 5c 4 t£ 4 )/( 
% 1) - (50 - 26 c 4 tl 3 ft? - 30 tl c 4 ft? 3 + 16 c 4 tl 2 t2 2 - 30 tl c 2 ft? -12 c 4 t2 3 tl 
+ 26 c 2 « *jg + 70 c 4 tl 2 ft? 2 + 10 t3 2 c 4 ft? 2 + 10 ft? ( 2 t2-8 ft? c 4 t2 3 - 62 c 4 tl 3 ft? 
injo , , n 4%3c 3 tlft£ 2 7%3c 3 tl 2 ft£ U%3<?tl 2 t3 

+ + ^i %i 

6°/c3c 3 tl t2t3 16 °/c3 tl c 3 ft? 2 4 %3 ft? c 3 tg 2 %Zt3 2 c?t2 . 4 4 
°/cl + %1 %1 " %1 + ° C U 

%3c« _ 5%3c 3 ft? 3 4°/c3ct# 2%3c 3 tl 3 _ 5%3cft? 2 2 
+ %1 %1 %1 %l %1 160 U 

+ 2 c 2 tl 2 + 46 c 4 tl 2 ft? ft? + 4 c 4 tl t2 2 ft? - 50 c 4 tl ft? ft? 2 + 10 c 2 ft? 2 + 22 c 4 tl 4 
+ 5c 4 ft£ 4 )/(%l)) + c 2 tl 2 - 2c 2 tl fcg + c 2 t^ 2 )^) 
%1 := 5c 2 ti? 2 + 2c 2 tl 2 + 5c 2 t3 2 - 9 - 2 tl c 2 ft? - 2c 2 tl ft? - 8 ft?c 2 ft? 
%2 := 50 - c° tl 2 tj? 4 + 12c 6 tl 5 ft? + 6c° tl 3 ft? 3 - 13c 6 tl 4 t2 2 + 12c 6 tl 5 ft? - c 6 tl 
+ 6c° tl 3 t3 3 - 13c 6 tl 4 ft? 2 - 16c 4 tl 3 ft? - 14c 6 tl 2 ft? 3 ft? + 34 c° tl 3 ft? 2 

- 34 c° tl 4 ft? t2 - t3 2 c 6 ft? 4 - 2 ft? 3 c° ft? 3 - ft? 4 c 6 ft? 2 + 34 c 6 tl 3 ft? 2 t2 

- 30 tl c 4 t3 3 + 8 c 4 tl 2 ft? 2 - 30 c 6 tl 2 ft? 2 ft? 2 - 20 tl c 2 ft? - 10 c 4 t£ 3 tl 

- 14 c° tl 2 ft? 3 t3 + 20 c 2 tl t2 + 10 c 6 tl t2 2 t3 3 + 2 c 6 tl ft? 4 ft? + 10 c G tl ft? 3 ft? 2 
+ 68 c 4 tl 2 ft? 2 + 2 tl c 6 ft? 4 ft? + 8 ft? 2 c 4 ft? 2 + 50 13 c 2 ft? - 10 13 c 4 ft? 3 - 4 c 6 tl 6 

- 56 c 4 tl 3 ft? + 10 13 3 c 4 ft? + 5 c 4 ft? 4 - 35 c 2 12 2 + 32 c 4 tl 2 t3t2 + 14 c 4 tl ft? 2 t3 

- 46 c 4 tl t2 t3 2 -15 c 2 ft? 2 + 18 c 4 tl 4 + 5 c 4 ft? 4 

%3 26 c 3 tl 2 t3 - 14 c 3 tl 2 ft? + 8 c 3 tl ft? 2 - 32 tl c 3 t3 2 + 2 ft? 2 c 3 t2 - 8 ft? c 3 t2 2 
- 8cft? + 10c 3 ft? 3 -2ctl + 10ct3 + 12c 3 tl ft? t3 - 4c 3 tl 3 - 6 V%2 
%4 := 26 c 3 tl 2 ft? - 14 c 3 tl 2 ft? + 8 c 3 tl ft? 2 - 32 tl c 3 ft? 2 + 2 13 2 c 3 ft? - 8 ft? c 3 t2 2 
-8ct2 + 10 c 3 1 ,? 3 - 2 c tl + 10 c 1 3 + 12 c 3 tl ft? ft? - 4 c 3 tl 3 + 6 V%2 

21. F- Value Of Target Position Report For Example #1 

> osmmts(tl J t2 J t3) [2] ; 

-=-{tl - t3) c((2 tl - ft? - t3) c + ((-(50 - - 26c 4 tl 3 ft? - 30 tl c 4 t, 

o /cl 

+ 16 c 4 tl 2 ft? 2 - 30 tl c 2 ft? - 12 c 4 ft? 3 tl + 26 c 2 tl ft? + 70 c 4 tl 2 t3 2 

+ 10 ft? 2 c 4 t2 2 + 10 ft? c 2 t2 - 8 « c 4 ft? 3 - %4 * + 16 %4 » ^ 

7c 1 Tel 

_ 4%4c 3 tl t2 2 7%4c 3 tl 2 ftg _ _ 6%4c 3 tl ft? ft? 

%1 + %1 %1 + 

4 °/ c 4 ft? c 3 tg 2 . 4 4 5%4cft? %4 ctl 5%4c 3 ft? 3 4°/ c 4ct^ 
+ %1 ~ + 0C td %1 + %1 %1 + %1 

2 %4 c 3 tl 3 

+ 0/1 -18c 2 t2 2 + 2c 2 tl 2 + 4d^tl 2 t3t2 + 4c i tlt2 2 t3 

/cl 
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or 



However, we have 

R ( x*-x 0 \ ( y*-y 0 \ 2 

(m. i -m 2 )(X( ! ;-xo) 2 +(l-^)(l/g-l/o) 2 

(m i -m 2 )(l-™ 2 ) 
(mi-m 2 )(x;-xo) 2 +(l-™ 2 )(i/S-!/o) 2 

((m 1 -m 2 ) 2 \ 

K - m 2 ) (a* - *o) 2 + (i - W - yof = { ^^R 

\ '"'1 / '"•l 

By Definitions 11 and 12, we have 

rrii 

Hence, the likelihood of the Target T being within the Error Likelihood Ellipse with constant ^ mi ~™ 2 ^ -R 
is calculated by the probabilities of M 3 (^.t^*!)- 
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51. (Mitigation Basis Algorithm 51 ) = 
mitig : — proc (xx-0 , yy-0 ) 

locals, x s.y s.d. d 1 ,d 2,g,s0,x0 s.yO s: 
globalx l.y l.x r.y r.x 1 ,y 1 : 

options ' Copyright 2003PQI Consulting All Rights Reserved': 

(ieschjrtion "OSMMTSuMitigationu Algorithm"; 

if (xx 0 , yy 0) — (x 1 , y 1 )iftenreturn —1,-1, —1. — lend if ; 

x s: — x I + s * (x r — x I): 

y s: - y I + s*(y r -y I): 

dl: = sqrt((((x 1 - x I) - s * (* r - x I)) © 2) + (((</ 1 - y I) - s * (y r - y I)) © 2)); 
d 2: - sqrt((((xx 0 - x I) - s * (x r - x I)) © 2) + (((yy 0 — y I) — s *(y r — y I)) ® 2)); 
d: - sqrt((x 1 - xx 0) © 2 + (y 1 - yy 0) © 2); 

g: — simplify ( 1 — 2 * (((((xx 0 — x I) — s * (x r — x I)) * ((1 — s) * (x r — x I)) 

+((W 0-y l)s*(y r-y l))*((l-s)*(y r-y 0))©2) / ((((xx 0-x l)-s*(x r-x l))®2+((yy 0 
-s * (y-r - yJ)) © 2) * (((1 - s) * (x.r - x.l)) © 2 + ((1 - s) * (y.r - yd)) © 2)) ) ) : 

sO : - fsolve(d-l © 2 + d-2 © 2 - 2 * d.l *d-2 *g = d, s. s — 0..1); 

if substring (convert (s 0 ., string). 1.. 3) <> "f so" then 

return sO,subs(s — sO .,xs). subs(s — sO.ys). subs(s — sO, d-1 + d-2); 

else return —1, minimize (g. s — 0..1), maximize(g. s — 0..1), —1; 

end if : 
This code is used in chunk 50. 

52. { Mitigation Extrema Algorithm 52 ) = 
xmit: — proc() 

local res , minres , maxres , u, v: 
globalx 0,y 0: 

options ' Copyright 2003PQI Consulting All Rights Reserved': 
rfe5cnptton"QSMMTSuMitigation u ExtreiiiauCalciilation": 

minres: — mitig (x-0 , y-0)[4]: 
maxres: — minres: 

for u from — 30to30 do for v from — 30to30 do res: — mitig(x 0 + u/10, y 0 + v/10)[4]; 

if res <> —lthen if res < minres then minres: — res end if : 

if res > maxres then maxres: — res end if : 

end if : 

end do : 

end do : 

return minres . maxres : 
end proc: 
This code is used in chunk 50. 

53. Sample Calculations And Results 

> (x_l,y_l):=(l,5); (x.r ,y_r) :=(2,2) ; (x_0,y_0) :=(4,3) ; 

> (x_l,y_l):=(3,6); 

xJ. yJ :— 1, 5 
x.r. y^r :— 2. 2 
x 0,y 0 := 4, 3 
x 1 , y 1 :— 3, 6 

> mitig(x_0-l,y_0+3) ; 

-1, -1, -1, -1 

> xmitO; 

3.493116110, 9.112507778 
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xx-0: 51. 
xO: 16. 
xOs: 51. 

xl : 16, 18, 23, 36, 39. 
x'2: 16, 18, 23, 36, 39. 
x3: 16, 18, 23, 36, 39. 
yd: 51. 
y.r; 51. 

: 51. 
y-0: 52. 
j/-/ : 51. 
yy-0: 51. 
2/0: 16. 
yOs: 51. 

16, 18, 23, 36, 39. 
y'2: 16, 18, 23, 36, 39. 
y3: 16, 18, 23, 36, 39. 



PQIC Documentation — 



OSMMTS Patent Application Technical Documentation 



11 



and 



L-2 



x — 



2(x 3 —x L ) 



J \ V3-Vl J 

for distinct Ki and k 2 , themselves functions of distinct Oi and 0 2 - 



±f2Cy/iK 2 +(f2cf+(x 2 3 -x^)+(y 2 3 -y 2 l ) 

2(3/3—3/1) 



Vi - 2/3 
2/1 ¥- 2/3 



do. 



III •- 



Figure 6: Third SDU Added 

However, since «u — sn (see the common side in Figure 6), then by (5), we have a link between Ki 
and k-2; namely 



or 



i(/ 1 c± V / 4 Kl + (/ 1 c) 2 ) 

(/i + / 2 ) c ± -/^i + (A c ) 2 ^ + ^ c ) 2 




x — 



±flC y /iK 1 + (f 1 cf+(x*-xT) 

2(x 2 —xi'j_ 

zb/ i c x /4/c 1 + (/ 1 c) 2 +(^-.rf) + (|/|-y;) 
2(1/2-1/1) 



y \V2-yi / 



and 



£2 : 



_ /2c((/ 1 +/2)c±V , 4/C i +(/ iC ) 2 )+(xj-X^) 
X ~ 2(X3-X 1 ) 

/2c((/i+/2)c±y , 4« 1 +(Ac) 2 ) + (.r 2 -.r 2 )+(i/ 2 -i/ 2 ) 
2 ('l/3—i/i) 



2/1 - 2/2 
2/1 ¥- 2/2 

2/1 = 2/3 
. 2/1 ^ 2/3 



(8a) 



(9a) 



v = -*(f^) 

Claim 3 If I. II, and III are not collinear, then Ly and L 2 are not parallel. 

Proof. Let /i^ be the slope of the line between SDU i and j. e.g., // 2 3 — ^ 2 I^ , or /x 2 3 — +00 for a 
vertical line. Since I, II, and III are not collinear, we have 

th2 ^ P23 j- I'rs 

with allowances for dealing with infinities, and Xi,X2,x$ are not all the same value. 

Suppose yi ^ y 2 and y± ^ y 3 . Then the slope of Ly is — ^ 2 ~^ ^ , and for L 2 it is — . If 



_ fx 2 - Xi \ _ _ fxs - Xi \ 

V 2/2 -2/1/ \ 2/3 -2/1/ 
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- 50 c 4 tl t'2 ft? 2 + 10 c 2 ft? 2 + 22 c 4 tl 4 + 5 c 4 t,? 4 ) /(7cl), -(50 - 26 c 4 £1 3 ft? 

- 30 tl c 4 ft? 3 + 16 c 4 tl 2 ft? 2 - 30 £1 c 2 ft? - 12 c 4 ft? 3 tl + 26 c 2 tl t2 + 70 c 4 tl 2 ft? 2 
+ 10 13 2 c 4 ft? 2 + 10 13 c 2 t2-S t3 c 4 ft? 3 - 62 c 4 £1 3 t3 + 10 ft? 3 c 4 t2 

_ 4 %3 c 3 tl t2 2 7 %3 c 3 tl 2 t2 _ 13 °/c3 c 3 £1 2 £3 _ 6 73 c 3 £1 t2 t3 
%1 + %1 7cl %1 

16 %3 tic? ft? 2 4%3ft?c 3 £g 2 %3t3 2 c 3 t2 . 4 4 73 cftl 
+ °/cl + 7cl %1 + ° C W + %1 

5 73 c 3 ft? 3 4%3_ctj? 2 7c3c 3 tl 3 5 73 eft? 1 „,,, 2 , 9 ,, r 2 
°/cl " + %1 + " %1 °/cl 18 C W +Z ° " 

+ 46c 4 £1 2 ft? ££ + 4c 4 £1 ft? 2 ft? - 50c 4 £1 ft? ft? 2 + 10c 2 ft? 2 + 22c 4 ft/ 4 + 5c 4 ft? 4 )/( 

%1)) + c 2 tl 2 - 2c 2 tl t2 + c 2 t2 2 )W) + 13 + J(« - t») c((-(50 

7c 1 

- 12 c 4 ft? 3 tl + 26 c 2 tl ft? + 70 c 4 tl 2 ft? 2 + 10 ft? 2 c 4 ft? 2 + 10 ft? c 2 ft? - 8 ft? c 4 ft? 3 

74 ft? 2 c 3 ft? 16 74 tl c 3 ft? 2 4 7c4 c 3 tl ft? 2 7 %4 c 3 tl 2 ft? m , 4 , „ 

- — 62 c 4 tl ft? 



7cl 7cl 71 71 

t™*™* + I0ft? 3 c 4 t2 + 4%4gc 3 t. 2 + , c4tg4 _ 574ct^ 
7c 1 7c 1 7c 1 

7c4ctl 5 7c4c 3 ft? 3 4%4cfc8 2%4c 3 tl 3 1D 2 . 2 „ . 2 . 2 

- 18c 2 t2 2 +2c 2 tl 2 



7cl 7cl 7cl 7cl 

+ 46c 4 tl 2 ft? ft? + 4c 4 tl ft? 2 ft? - 50c 4 tl £2 ft? 2 + 10c 2 ft? 2 + 22c 4 tl 4 + 5c 4 ft? 4 )/( 

7c 1), -(50 - 26 c 4 tl 3 ft? - 30 tl c 4 ft? 3 + 16 c 4 £1 2 £2 2 - 30 tl c 2 ft? - 12 c 4 ft? 3 tl 
+ 26 c 2 tl ft? + 70 c 4 tl 2 ft? 2 + 10 ft? 2 c 4 ft? 2 + 10 ft? c 2 £2 - 8 ft? c 4 £i? 3 - 62 c 4 tl 3 ft? 

0 , 4 „ 473c 3 £1 ft? 2 773c 3 £1 2 ft? 137c3c 3 £l 2 ft? 
+ 10t?C ^- 71 + 71 - 71 
_ 6 73 c 3 tl ft? ft? 16 73 tl c 3 ft? 2 4 73 ft? c 3 £2 2 _ 73 ft? 2 c 3 ft? . 4 4 
71 + 71 + 71 71 + ° C 

%Zctl 5 73c 3 ft? 3 4 73 eft? 273c 3 tl 3 5%3ctg 2 2 
+ 71 71 + 71 + 71 71 8C U 

+ 2 c 2 tl 2 + 46 c 4 tl 2 ft? ft? + 4 c 4 tl ft? 2 ft? - 50 c 4 tl £2 £3 2 + 10 c 2 ft? 2 + 22 c 4 tl 4 
+ 5 c 4 t^ 4 )/(7l)) + c 2 tl 2 - 2 c 2 tl ft? + c 2 £^ 2 )(V2) 

71 := 5 c 2 ft? 2 + 2 c 2 tl 2 + 5 c 2 ft? 2 - 9 - 2 tl c 2 ft? - 2 c 2 £1 ft? - 8 ft? c 2 ft? 

72 := 50 - c G tl 2 ft? 4 + 12 c G tl 5 ft? + 6 c° tl 3 ft? 3 - 13c 6 tl 4 ft? 2 + 12 c G tl 5 ft? - c 6 tl 2 ft? 4 
+ 6c° tl 3 ft? 3 - 13c 6 tl 4 ft? 2 - 16c 4 tl 3 t2 - 14c 6 tl' 2 ft? 3 ft? + 34c 6 tl 3 t3 ft? 2 

- 34 c 6 £1 4 ft? ft? -t3 2 c°t2 i -2 ft? 3 c 6 ft? 3 - ft? 4 c 6 ft? 2 + 34 c° £1 3 ft? 2 £2 

- 30 tl c 4 ft? 3 + 8 c 4 £1 2 ft? 2 - 30 c 6 tl 2 t3 2 t2 2 - 20 £1 c 2 ft? - 10 c 4 ft? 3 tl 

- 14 c 6 £1 2 ft? 3 ft? + 20 c 2 £1 ft? + 10 c 6 tl ft? 2 ft? 3 + 2 c 6 £1 ft? 4 ft? + 10 c 6 tl t£ 3 ft? 2 
+ 68 c 4 £1 2 ft? 2 + 2 tl c 6 ft? 4 ft? + 8 ft? 2 c 4 ft? 2 + 50 ft? c 2 ft? - 10 £5 c 4 ft? 3 - 4 c 6 tl 6 

- 56 c 4 £1 3 ft? + 10 ft? 3 c 4 ft? + 5 c 4 ft? 4 - 35 c 2 ft? 2 + 32 c 4 tl 2 ft? ft? + 14 c 4 tl t£ 2 ft? 

- 46 c 4 £1 ft? ft? 2 -15 c 2 ft? 2 + 18 c 4 £1 4 + 5 c 4 ft? 4 
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31. Error Likelihood Ellipse Calculation 

The Error Likelihood Ellipse is the totality of Target Positions Reports for a particular set of three SDU ; s 
produced by varying the observed arrival time data according to an assigned error distribution. The resulting 
(x, y) position that has maximum likelihood within this ellipse is judged to be the most likely "real" position 
of the target. The same algorithms used to calculate an individual Target Position Report shall be used to 
calculate the Error Likelihood Ellipse. 

32. Error Likelihood Ellipse Standard Methodology 

The PQIC OSMMTS Error Likelihood Ellipse Standard Methodology is to assume the arrival time data 
varies according to a normal distribution of mean zero and standard deviation a, where a is either estimated 
from calibration data or the result of a Bayesian estimation process. This means the observed arrival time 
data is always accepted as the center of the error distribution. 

The algorithms found below implement such an error distribution. However, such algorithms may be 
amended at any time to more accurately reflect the true error distribution in any particular PQIC OSMMTS 
System implementation. 

33. Error Likelihood Ellipse MAPLE Algorithm 

Given the position information of the SDU's that receive the arrival time data, the following MAPLE 
algorithm calculates the Error Likelihood Ellipse based on the normal error distribution affected arrival 
times. These algorithms call the Target Position Report algorithms by reference, so that the Error Likelihood 
Ellipse algorithms are independent of the particular methods found in the Target Position Report code. 

The M 0 distribution and the Target Position Report algorithms are needed to calculate the Error 
Likelihood Ellipse. All other M distributions are only needed when quantifying the bounds of the error 
ellipse, or for making multivariate calculations on multiple Target Position Reports and Error Likelihood 
Ellipses. 

34. { ELE Calculation 34 ) = 
MOGEN: - proc {tl,t2, tS , a, b, r) 
local ans ,el,e2, e3 , jl , j2 , j3 ■ idx , tx : 

options ' Copyright'2003PQI Consulting All Rights Reserved': 
rfescripiwn "OSMMTS u MOuGeneration u Cal dilation": 
idx: — 0; 
with (stats): 

for el fromatob do jl : — statevalf [cdf . normald[0.0.1]](el /r) — statevalf [cdf , normald[0,0,l]]((el — l)/r): 
for e2 fromatob do j2: — statevalf [cdf , normald[0, 0.1]] (el? /r) — statevalf [cdf , normald [0,0.1]] ((e2 — l)/r): 
for e3 fromatob do j3: — statevalf [cdf ,normald[0,0.l]](e3 /r)— statevalf [cdf , normald[0,0.1]]((e3 — l)/r): 
tx: - evalf(osmmts(tl + el jr. t2 + e2/r, t3 + e3/r)[],b): 
if (type(tx[S\, undefined)) then 
else 

if (abs(tx[4]) > 1) then 
else 

idx: — idx + 1; 

ans[idx]: = (tx[l],tx[2],jl *j2 *j3): 

end if : 

end if : 

end do : 

end do : 

end do : 

return ans , idx : 

end proc: 
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54. Optimal Placement Of Detection Units 



/ ♦ T II 



! V ! 
/// * > IV 



Figure 11: Poor Placement Of SDU Units 

The optimal placement of SDU's depends on the OSMMTS goals for the implementation under consid- 
eration. In general, a maximum likelihood approach to SDU placement works well for many applications. 
Its goals address the probability of an accurate TPR calculation with the first qualifying set of SDU's within 
a particular SDU arrangement. The probability of such placement is maximized among all convex polygons 
that encloses all possible target locations, or those target locations that are of interest in the implementation, 
or where the addition of any more SDU's would not enclose a significantly larger area than is covered by 
smaller set of SDU's. Such an approach uses the following rules for calculating a set of exact SDU positions 
for optimal use. 

Optimized SDU Placement Rules 

1. SDU's should be placed where the connected units form a convex hull (see Figure 12 for an example of 
a convex hull, and Figure 11 for an example of a non-convex hull). 

2. The convex hull should cover as much of the possible target locations as possible. 

3. SDU's are placed so that no three units are colhnear. 

4. SDU's are placed so that no two units have the same avralue or y-value, as viewed in a coordinate grid. 

5. SDU's should be placed to maximize the likelihood that di > \ fi\c. for all i. 

No consideration is made in these rules for obstructions, or other exceptions in a clear area of target 
detection. Under the OSMMTS, such considerations are handled by mitigations. Furthermore, if the set of 
possible target locations, or those target locations that are of interest in the implementation, contain "gaps," 
regardless of size and shape, such gaps shall not be considered in the optimal placement of detection units, 
since their consideration may be completely addressed by the methods found in the containment policies. 
Specialized implementation of the OSMMTS may introduce special SDU placement considerations, which 
will be documented under separate cover in the particular implementation documentation. 

Figure 13 shows an arrangement of five SDU's that satisfy the guidelines embedded in the placement 
rules. Since the position of the origin is arbitrary, the coordinate center of the region may be translated to 
any local coordinate system within any particular OSMMTS implementation. 
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( ELE Example #1 SDU Coordinates 36 ) Used in chunk 12. 

( ELE Example #2 SDU Coordinates 39 } Used in chunk 12. 
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The Optimal Surface Mitigated Multiple Targeting System 
Patent Application Technical Documentation 

Timothy Hall 

Abstract 

The purpose of the PQIC Optimal Surface Mitigated Multiple Targeting System (PQIC 
OSMMTS) is to define the equipment and processing necessary to produce, in real time, 
an error-bounded, self-monitoring and self-adjusting, likelihood-based Target Position 
Report for arbitrarily many self-identifying targets in a two-dimensional grid. Each 
target sends identifying information to an array of sensors strategically placed in its 
vicinity to maximize the likelihood that the system will calculate a position report as 
accurately and precisely as possible. 
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Figure 7 : Intersection of Foci Locates T 



then either 

2/2 - 2/1 2/3 - 2/1 



or 



or 



x 2 —xi x 3 — Xi 

2/3 - 2/i , 

+oc — . when x± — x 2 

x 3 -x±' 

2/2 - 2/1 



+oc. when x± — x 3 

X 2 - Xi 

E x*-x\ ~ X3-X1 ; ^en /ii2 — M13; a contradiction. 

If +oc — v ?~ v } . then x\ — x 3 . since y± / w 3 . This is also a contradiction, since x± — x 2 under this 
situation. 

E ^I^i ~ + oc - then a?i = » 2 » since ?/i ^ jfe, yet another contradiction, since xi — x 3 under this 
situation. 

If either yi — y 2 and y 2 ^ y 3 , or y 2 — y 3 and 1/1 ^ j/ 25 then the slope of L\ would be +oc, and for L 2 it 
would be — ( ^Zyl ^-, or the other way around. In either case, these slopes could not be equal in any sense. 

Finally, if both y± — y 2 and y± — y 3 , then I, II, and III would be collinear, since they would all have 
the same y-value: this is also a contradiction. ■ 

By Claim 3, let 9 0 be the value of 9 for the point of intersection between L\ and L 2 . When yi 7^ y 2 
and yx ± y 3 . let 

' X 3 - Xi. 



__( X-2 ~ Xl \ 

\y2-y1) 



nil — — — and m 2 — — 



2/3 - 2/1 



with 



for kq corresponding to 9 0 . 
Then 



b _ ±ficy/4Kp + (Ac) 2 + (4 - 4) + (2/f - yj) 
2(j/2-2/i) 

^ _ ±f 2 cyj4Ko + {f2cf + (4 - 4) + {yi - y'l) 

2 (y 3 - j/i) 



Li : y — mix + bi 
L 2 :y — m 2 x + b 2 
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%3 := 26 c 3 tl 2 t3 - 14 c 3 tl 2 t2+8 c 3 tl t'2 2 - 32 tl c 3 t^ 2 + 2 t,? 2 c 3 ft? - 8 t3 c 3 ti? 2 

- 8cti? + 10c 3 t,? 3 -2ctl + I0ct3 + 12c 3 tl t2 W - 4c 3 tl 3 - 6 V%2 

%4 := 26 c 3 tl 2 i3 - 14 c 3 tl 2 t2 + 8 c 3 « t? 2 - 32 tl c 3 t3 2 + 2 tf 2 c 3 t2 - 8 *5 c 3 fc? 2 

- 8 c ftg + 10 c 3 tS 3 - 2 c « + 10 c + 12 c 3 tl t2 tS - 4 c 3 tl 3 + 6 V%2 

22. Target Position Report Example #2 

23. { TPR Example #2 SDU Coordinates 23 ) = 
{xl,yl): = (3,2); (x2,y2): = (2,4); (atf,y5): - (1,3); 

This code is used in chunk 12. 



24. Target Position Report MAPLE Calculation For Example #2 

> c:=l; (xl,yl):=(3,2);(x2 J y2):=(2 J 4);(x3,y3):=(l J 3); 

c:= 1 
ari, yl := 3, 2 
^, y£ := 2, 4 
atf, := I, 3 

> osmmts(tl J t2 J t3) [1] ; 

25. A- Value Of Target Position Report For Example #2 

-\{tl - t3) c((2 tl-t2- t3) c + (4min(-(50 + 26c 2 tl t£ - 8 W c 4 t2 3 - 30 tl c 2 t3 

O 

r n,o 2,0 4 6%Zc 3 tlt2t3 4%3t3c 3 t2 2 %3t3 2 c 3 t2 

+ 10 t,? C 2 t£ + 10 t,? 3 C 4 fcg ^ h ft- ft- 

% 1 7cl 7cl 

T /:)<■'■ il 2 i ! 4%3c^tlt2 2 W%3tlc?t3 2 13 %3 c 3 tl 2 t,? 

+ 7^r- 7^r- V 



%1 %1 %1 %l 

tl 3 t 
o%3ct3 



. 62c 4 tl 3 2% 1 C ! ~ ~ 5% tt ^ - » 3 * + ^ 

7C 1 Tel 7c 1 7c 1 



+ 10 t3 2 c 4 ti? 2 - 30 tl c 4 t,? 3 + 70 c 4 tl 2 t3 2 - 12 c 4 t2 3 tl 

7c 1 

+ 16 c 4 tl 2 t2 2 + o c 4 t3 4 - 18 c 2 ti? 2 + 2 c 2 tl 2 + 46 c 4 tl 2 t3t2 + 4 c 4 tl t£ 2 t5 
- 50 c 4 tl t2 t3 2 + 10 c 2 t3 2 + 22 c 4 tl 4 + 5 c 4 t2 A ) /(4 %1), -(50 + 26 c 2 tl fc? 

4%4t,?c 3 ^ 2 



- 8 13 c 4 ti? 3 - 30 tl c 2 £? + 10 13 c 2 i? + 10 t3 3 c 4 fc? + 



%1 



c /c4 tg 2 c 3 t2 7 %4 ( s tl 2 t2 4%4c 3 tlt2 2 16 %4 tl c 3 t,? 2 

c /cl + %1 %1 + 7cl 

_ iVM^tl _ ,%4^1t2t3 _ 62c4 tl3 t3 _ 26c4 tf 1Q ^ c4 ^ 

7c 1 Vol 
5 0 /c4c 3 tff 3 4%4ct2_ _o%4ct3_ %4ctl 2%4c 3 « 3 4 3 

°/cl ^ + %1 %1 + °/cl + °/cl 

+ 70 c 4 tl 2 tg 2 -12 c 4 ti? 3 tl + 16 c 4 tl 2 t2 2 + 5 c 4 t,? 4 - 18 c 2 t2 2 + 2 c 2 tl 2 
+ 46 c 4 tl 2 t3 t2 + 4 c 4 tl t2 2 tS - 50 c 4 tl ft? tg 2 + 10 c 2 t,? 2 + 22 c 4 tl 4 + 5 c 4 ti? 4 )/( 

4 %1)) + c 2 tl 2 - 2 c 2 tl tl? + c 2 t2 2 fW) + ^ + \i tl ~ t2 ) c(4min(-(50 
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J: — proc(dt , cnt , msh) 

local idx , inx , iny , new , smm , minx , maxx , miny . maxy ,g,gx,gy,M: 
options ' Copyright'2003PQI Consulting All Rights Reserved': 
description "OSMMTS u MO u Processing" ; 
minx: — 0; maxx: — 0; 
miny: — 0; maxy: — 0; 
9- -0; 

for idx fromltocnt do new: — dt[idx]:. 
inx: — floor (msh * new[l]); 
iny: — floor (msh * new [2]): 
if (assigned (M [inx, iny])) then 
M[inx,iny]: — M[inx, iny] + new [3]; 
else M[inx,iny]: — new [3]; 
end if : 

if (inx > maxx) then 
maxx: — inx: 
end if : 

if (inx < minx) then 
minx: — inx: 
end if ; 

if (iny > maxy) then 
maxy: — inx: 
end if : 

if (iny < miny) then 
miny: — inx: 
end if : 
end do : 

for idx from minx to maxx do for iny from miny to maxy do 

if (assigned (M [idx, iny])) then 
if (M [idx, iny] > g) then 
g: — M[idx,iny]: 
gx: — idx: 
gy: - iny: 
end if : 

else M[idx,iny]: — 0: 
end if : 
end do : 
end do : 

return eval (M), minx , maxx , miny , maxy .gx.gy.g: 
endproc: 
This code is used in chunk 12. 

35. Error Likelihood Ellipse Example #1 

36. (ELE Example #1 SDU Coordinates 36) = 
(xl,yl): = (2,3); (x2,y2): - (4,2); (x3,y3): - (3,1); 

This code is used in chunk 12. 

37. Error Likelihood Ellipse MAPLE Calculation For Example #1 

> with(stats) : 

> mesh: =35; 

> TPR : = (evalf (osmmts (16/10,2, 175/100) [1 . . 2] , 10) ) ; 

TPR := 2.547503694, 2.097308261 
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Figure 12: Better Placement Of SDU Units 
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Figure 13: Best Placement Of SDU Units 
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Claim 4 When I. II. and III are not collinear. and yi ^ y 2 and 1/1 ^- y ?> , the point of intersection between 
Ll and,L 2 is (* 0 ,2/o) = (^.ffn (^) + &i). 
Proof. We have 

mixo + bi — m 2 x 0 + f> 2 



or 



and 



Xq — 



b-2-h 
mi — m 2 



\ Li - mxxo + 61 

= ™i(^)+&i 

= (m 1 -m 2 + m 2 )( 7 ^) + & 1 

= 62 _ 6l + 6l+m2 (^_) 

- m 2 (^)+t>2 



l£ 2 



Claim 5 When I, II, and III are not collinear. and y± — y 2 and yi 7^ y 3 . the point of intersection between 
Li and L 2 is 

( 



(#0,2/0) 



2(x 2 —xi) 



2(3/3-2/1) 



V 

Proof. We have 



he ({fi + h)c± (Ac) 2 ) + (4 - x\) + (yi - yf) ^ 

- (±hcj4Ko + (fic) 2 + (xl-xl)^ (g : ) j 



Li : x — 



±flC^4Kp + (he) 



A 2 4- /V? 



K — X 



1) 



and 



L 2 :y = -x 



x 3 - Xi 
2/3 - 2/1 



+ 



Hence 
2/0 1 i2 - 



2 (# 2 - #1) 

he ((A + A) c + ^4^ + (Ac) 2 ) 
2(2/3-2/1) 



(2/1 - 2/ 2 ) 



+ 



_ ( ±/ 1 c v / 4«o+(/ic) 2 +(xl-x?) 
2(x 2 -xi) 



/ ±/icV , 4/t 0 +(/ic_) 2 +(a;|-x?) 
I 2(3/3-3/1) 



/2c((/ 1 +/ 2 )c± A /4to+(/ic) 2 )+(xl-x 2 )+(y|-;/ 2 ) 
2(3/3-3/1) 

/2 C ((/i+/2)c±v'4 re "o + (/ lC ) 2 ) + (x 2 -x 2 ) + (,/ 2 -;/ 2 ) 

2(3/3-3/1) 

/2c((/l + /2)c±V'4«0 + (/l C ) 2 )+(x 2 -X 2 ) + (,/ 2 -,/ 2 ) 

2(3/3-2/1) 



X3- 
3/3" 



2(3/3-3/l) 



-XI ^ 

-3/1 / 

\ Z' X 3 -Xj\ 

/ \X2-X1 / 



f he ({f i + h) c ± Jako + (Ac) 2 ) + (a§ - s?) + (2/1 - 2/ 2 ) X 

^ - (iAc^Ko + (Ac) 2 + (4 - 4)) (:iz:i 



Claim 6 When I, II, and III are not collinear, and yi ^ y 2 and yi — 2/3 , the point of intersection between 
Li and L 2 is 

( /2C 2 (/i + / 2 )±V'4ko+(/ic) 2 +(x 2 -x 2 ) \ 



(£0,2/0) 



2(2/2-2/1) 



/ 2 c 2 (/i + / 2 )±\/4/vo+(/ic) 2 +(x 2 -.i- 2 ) 
2(x 3 -x i) ' 

±Acy / 4'Co + (Ac) 2 + (*| - xf) + (2/I - J/ 2 ) 

/ 2 c ((a + A) c ± v / 4k 0 + (Ac) 2 ) + (*§ - *?)) (:: 2 ;;: ) 
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+ 26 c 2 tl ft? - 8 t3 c 4 ft? 3 - 30 tl c 2 ft? + 10 ft? c 2 ft? + 10 ft? 3 c 4 ft? 
6 %3 c 3 tl ft? ft? 4%3t3c 3 t2 2 _ 73 ft? 2 c 3 ft? 7 %3 c 3 tl 2 ft 



7cl %1 71 7cl 

4 %3 c 3 tl tg 2 , 16 7c3 « c 3 ft? 2 13 7c3 c 3 tl 2 ft? „„ 4 3 , „ , %3 c tl 



+ ^ ^ 62 c A trt3 + 



%1 7cl 7cl %1 

2%3c 3 tl 3 5 73 c 3 ft? 3 4, r j,„ 4 73 eft? 573ct,? 

+ c /cl %1 ° C W + %1 %1 

+ 10 13 2 c 4 t2 2 - 30 tl c 4 ft? 3 + 70 c 4 tl 2 ft? 2 - 12 c 4 ft? 3 tl + 16 c 4 tl 2 ft? 2 + 5 c 4 ft? 4 

- 18 c 2 t£ 2 + 2 c 2 tl 2 + 46 c 4 tl 2 ft? ft? + 4 c 4 tl ft? 2 ft? - 50 c 4 tl tg ft? 2 + 10 c 2 ft? 2 
+ 22 c 4 tl 4 + 5 c 4 ft? 4 ) /(4 7c 1), -(50 + 26 c 2 tl ft? - 8 ft? c 4 ft? 3 - 30 tl c 2 ft? 

+ io ft?c 2 * + ioft? 3 c 4 1 2 + ******* _ + r%^ ftg 

/cl /cl /cl 

47c4c 3 tl tg 2 16 74 tl c 3 ft? 2 13 7c4 c 3 tl 2 ft? 6 74 c 3 tl ft? ft? 
71 + 71 71 71 

- 62c 4 tl 3 ft? - 26c 4 tl 3 t2 + 10 ft? 2 c 4 ft? 2 - + 4%4c ^ 5%4c ^ 



71 1 71 71 
V> [€il 2 " lj3 - 30 tl c 4 13 3 + 70c 4 tl 2 t3 2 - 12c 4 t2 3 tl + 16c 4 tl 2 ft? 2 



71 71 

+ 5 c 4 ft? 4 - 18 c 2 ft? 2 + 2 c 2 tl 2 + 46 c 4 tl 2 ft? ft? + 4 c 4 tl ft? 2 ft? - 50 c 4 tl t2 ft? 2 
+ 10c 2 ft? 2 + 22c 4 tl 4 + 5c 4 ft? 4 )/(47l)) + <? tl 2 - 2c 2 tl ft? + c 2 ft? 2 )^ 

71 := 2c 2 tl 2 + 5c 2 ft? 2 + 5c 2 ft? 2 -9 - 8ft? c 2 ft? - 2tl c 2 ft? - 2c 2 tl t2 

72 := 50 + 20 c 2 tl ft? - 10 ft? c 4 ft? 3 - 20 tl c 2 ft? + 50 ft? c 2 ft? + 10 ft? 3 c 4 ft? - 56 c 4 tl 3 ft? 

- 16 c 4 tl 3 t2 + 8 ft? 2 c 4 t2 2 - 30 tl c 4 ft? 3 + 68 c 4 tl 2 ft? 2 - 10 c 4 ti? 3 tl + 8 c 4 tl 2 t2 2 
+ 5 c 4 ft? 4 - 35 c 2 ft? 2 + 32 c 4 tl 2 ft? ft? + 14 c 4 tl ft? 2 ft? - 46 c 4 tl ft? ft? 2 -15 c 2 ft? 2 
+ 18 c 4 tl 4 + 5 c 4 ti? 4 - 4 c G tl c - 34 c G tl 4 ft? t2 + 34 c 6 tl 3 ft? ft? 2 - 14 c c tl 2 ft? 3 ft? 

- 30 c G tl 2 ft? 2 ti? 2 + 34 c G tl 3 ft? 2 ft? - 14 c c tl 2 ft? 3 ft? + 6 c 6 tl 3 ft? 3 + 12 c c tl 5 ft? 

- c G tl 2 ft? 4 - 13 c G tl 4 ft? 2 + 6 c c tl 3 ft? 3 + 12 c 6 tl 5 ft? - c 6 tl 2 ft? 4 - ft? 4 c G ft? 2 

- 2 ft? 3 c G ft? 3 - ft? 2 c G ft? 4 - 13 c G tl 4 ft? 2 + 10 c 6 tl ft? 3 ft? 2 + 2 c G tl t2 4 ft? 
+ 10 c G tl ft? 2 ft? 3 + 2 tl c G ft? 4 ft? 

73 := 26 c 3 tl 2 ft? - 14 c 3 tl 2 ft? + 8 c 3 tl ft? 2 - 32 tl c 3 ft? 2 + 2 ft? 2 c 3 ft? - 8 ft? c 3 ft? 2 

- 4c 3 tl 3 + 10c 3 ft? 3 + lOcft? - 8ct2 -2ctl + 12c 3 tl ft? ft? + 6 vT~2 

74 := 26 c 3 tl 2 ft? - 14 c 3 tl 2 ft? + 8 c 3 tl ft? 2 - 32 tl c 3 tf 2 + 2 ft? 2 c 3 ft? - 8 ft? c 3 ft? 2 

- 4 c 3 tl 3 + 10 c 3 ft? 3 + 10 c ft? - 8 c ft? - 2 c ft/ + 12 c 3 tl ft? t,? - 6 V%2 

26. F- Value Of Target Position Report For Example #2 

> osinmts(tl J t2,t3) [2] ; 

- g (tl - ft?) c((2 tl - ft? - ft?) c + (4min(-(50 + 26c 2 tl ft? - 8 ft? c 4 ft? 3 - 30 tl <? t3 

+ lot.? c 2 ft? + lot.? 3 c 4 1* - 6%3c i^ + 4%3 ;y^ 2 - 

/cl /cl /cl 

773c 3 ftl 2 ft? 4%3c?tlt2 2 16 73 tl c 3 ft? 2 13 73 c 3 tl 2 ft? 
+ 71 71 + 71 71 

- 62c 4 tl 3 ft? + **£L + 2%3 / " 3 - 5%3 / - 26c 4 tl 3 ft? + 

/cl /cl /cl /cl 
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> RES : =M0GEN (16/10,2, 175/100 , -12 , 12 ,mesh) : 

Warning, these names have been redefined: anova, describe, fit, 
importdata, random, statevalf, statplots, transform 

> MZ:=J(RES[l],RES[2],mesh): 

38. Error Likelihood Ellipse Example #2 

39. { ELE Example #2 SDU Coordinates 39 ) = 
{xl,yl): = (3,2); (x2,y2): = (2,4); (x3,y3): = (1,3); 

This code is used in cliunk 12. 

40. Error Likelihood Ellipse MAPLE Calculation For Example #2 

> with(stats) : 

> mesh: =35; 

> TPR : = (evalf (osmmts (16/10,2, 175/100) [1 . . 2] , 10) ) ; 
TPR := 2.097308261, 2.547503694 

> RES : =M0GEN (16/10,2, 175/100 ,-12,12 ,mesh) : 

Warning, these names have been redefined: anova, describe, fit, 
importdata, random, statevalf, statplots, transform 

> MZ:=J(RES[l],RES[2],mesh): 
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55. The Convex Hull Principle 



Since the area of any segmented convex hull is maximized when the distances between the vertices of the 
segments are equal 14 , the "best" convex hull, i.e., that which follows the placement rules, is based on a circle 

14 Since all convex polyhcdra may be subdivided into triangles, it suffices to prove the following claim to make this assertion. 
Claim 23 Given a fixed perimeter length, the triangle with maximum area is equilateral. 

Proof. Let a > 0, b > 0, and c > 0 be the length, of the sides of a triangle. Then the area of such a triangle is given by 

A — \/« (« — a) (s — b) (s — c) 

where s — | (a + b + c). It remains to show that A achieves its maximum when a — b — c, given a fixed value of a + b + c. 



Consider the Lagrange Multiplier Function K: 



K (a,b,c) — A — \ (a + b + c) 



Then 



where 



dK 

da 



dA 
da 



A: 



Bi 

4/1 



{-a + b + c)(a-b + c)(a + b- c) 
- (a + b + c) (a - b + c) (a + b - c) 
+ (a + b + c) (-a + b + c)(a + b- c) 
+ (a + b + c) (-a + b + c) (a - b + c) 



Similarly, we have 



where 



dK 
db 



dA 
db 



B 2 
4A 



B 2 ■■ 



Setting 



we have 



(-a + b + c)(a-b + c)(a + b-c) 
+ {a + b + c)\a-b + c)\a + b- c) 
-(a + b + c) (-a + b + c)(a + b-c) 
+ (a + b + c) (-a + b + c) (a-b + c) 



dK _ dK 
da db 



• (a + b + c) (a - b + c) (a + & • 



4/1 



B 2 

4/1 



which means 
Furthermore, 

gives 

or 

which means 
In fact, 

where 



Then a — b — c gives 



+ (a + b + c) (-a + b + c)(a + b 

a — b-\-c— —a + b + c 
a = & 



+ (a + & + c) (a - & + c) (a + b - c) 
-(a + b + c) (-a + b + c) (a + b - c) 



: 0 : 



dK 
dc 



f - (a + & + c) (a - 6 + c) (a + & - c) If + (a 
\ +(a + b + c){-a + b + c)(a-b + c) J \ -(o- 



+ b + c) (a - b + c) (a + b - c) 1 
+ & + c)(-a + <i + c)(a-<i + c) J 



o + 6- 



da 2 



-a + b + c 



4ABu - B l 
16/1 2 



-2(a-b + c)(a + b- c) 
+4a (-a + b + c) 
+2 (a + b + c) (-3a + b + c) 



and 



-4<r 



3 \ 1 /; A ^3 2 

a — a = V3a' = a 
2/4 4 
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1. Mission Statement 

The purpose of the PQIC™ Optimal Surface Mitigated Multiple Targeting System (PQIC OS- 
MMTS 1 ) is to define the equipment and processing necessary to produce, in real time, an error-bounded, self- 
monitoring and self-adjusting, likelihood-based Target Position Report for arbitrarily many self-identifying 
targets in a two-dimensional grid 2 . Each target sends identifying information to an array of sensors strate- 
gically placed in its vicinity to maximize the likelihood that the system will calculate a position report as 
accurately and precisely as possible. The OSMMTS uses analytical and ad-hoc mitigation and optimization 
techniques to reduce the error bounds on the target report to a practical minimum. This memorandum 
provides the baseline technical documentation of the analytical methods, construct guidelines, quantifi- 
cation methods, mitigation and optimization techniques, and programming details for implementing the 
PQIC OSMMTS. Additional supporting methods, documented separately, that supplement, complement, 
and refine the analytical methods found herein shall be issued in an ongoing series through PQIC technical 
documentation. 



0 Process Quality Improvement Consulting, PQIC, Optimal Surface Mitigated Multiple Targeting System, and OSMMTS 
arc trademarks of PQI Consulting, P. O. Box 425616, Cambridge, MA 02142-0012 USA, info@pqic.com." 

1 While the formal name of the system is PQIC OSMMTS for copyright, trademark, service mark, and patent purposes, 
the name most often used in this, and subsequent technical, marketing, and promotional documentation shall be the OSMMTS 
System, or simply the OSMMTS. 

2 A separate, yet related, targeting system for a three-dimensional grid will be covered in a separate set of PQIC documen- 
tation. 



14 



PQIC Documentation 



— OSMMTS Patent Application Technical Documentation 



Proof. We have 

Li : y - -x 



and 



Hence 



X 2 -Xi 

y-2 - 2/1 



+ 



i/iC^Ko + (Ac) 2 + (*! - xf) + (yl - y{) 
2(i/2 -1/1) 



L2 : x — 



he (ifi + h) c±^Ako+ (Ac) 2 ) + (4 - xf) 
2(x 3 -xi) 



--( 



2(3/2-3/1) 

M f ^ 

2(x 3 -xi) / \ 3/2-3/1 / 



_ / /2c((/i+/2)c±y , 4«o+(/ 1 c) 2 )+(xj-x^ > 



+ 



/2c((/ 1 +/ 2 )c± v / 4^+(7^F)+(x|-x?) 
2(3/2-3/1) 



/ X2— xi 
^x 3 — xi 



±h ct/4k 0 + ( /i c) 2 + (xj -x| ) + ( ;/ 2 -j/ 2 ) 
2('i/2— !/i ) 

±/ 1 cv'4/to+(/ic) 2 +(^-a- 2 )+('y 2 -j/ 2 )" 
2(lte-yi) 



2(3/2-3/1) 



V 



i/ic^Ko + (Ac) 2 + (a§ - #?) + (y2 - J/2) 
- ( Ac ('(A + A)c± y/bio + (Ac) 2 ) + (*§ - xf)) (%Z%l) 



The fourth possibility consistent with Claims 4-6, namely 1/1 — 1/2 and 1/1 — 2/3, has been excluded by 
the assumption that I, II, and III are not collinear. 

Note how in all Claims 4-6, (xq, y 0 ) is a function of constants and the Sy (via the bi), which are in turn 
functions of 9. This means 9 0 may be expressed implicitly through the point (xo,yo). 

If 

Vl = (#1,2/1) - (#0,2/0) 

and 

V 2 - (#2,2/2) - (#0,2/0) 

then |VI| 2 — sn and |I'2| 2 — «2i, so that 

VI • V 2 - |VI| 2 \V2\2 COS^o - »11«21 cosflo 

This proves the following claim. 

Claim 7 <?o = Arc cos f 

u \s11s21 

The quadrant of # 0 must be considered when calculating the inverse cosine in Claim 7. Once established, 
we then have 

Sll (A) S2I (A) COS # 0 - (#1 - #0 (0 O )) (#2 - #0 (A))) + (j/l ~ V0 (A))) (ife ~ M> (A)) (10) 

However, from (2), 

an (^o ) S21 (A)) - tk> 

and. from (3). 

COS *°- 1 2^" 

so 

an (A>) ^21 (A>) cos A) - «o ^ - ^ ~^ lC)2 ^ (11) 

or 

Ko = <*i-(A c ) + _ ^ 0 (0 O )) _ ^ (0 O )) + ( yi _ j/ 0 (0 O )) ( y2 _ j/ 0 (0 O )) 

by (10). 

Solving for k 0 . and substituting tliis value into Claim 4, 5, or 6, depending on circumstances, determines 
the position of Target T. In fact, 



* = *i - 7 = *i - ^ (Aciv^o + C^) 2 ) 
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- ° y g ctg + 10 ft? 2 c 4 ft? 2 - 30 tl c 4 ft? 3 + 70 c 4 tl 2 ft? 2 - 12 c 4 tl> 3 tl 

7c 1 

+ 16 c 4 tl 2 ft? 2 + 5 c 4 ft? 4 - 18 c 2 ft? 2 + 2 c 2 tl 2 + 46 c 4 tl 2 t3 t2+4 c 4 tl t2 2 t3 

- 50 c 4 tl ft? ft? 2 + 10 c 2 ft? 2 + 22 c 4 « 4 + 5 c 4 ft? 4 ) /(4 %1), -(50 + 26 c 2 tl t2 

- 8 ft? c 4 ft? 3 - 30 « c 2 t3 + 10 ft? c 2 t2 + 10 ft? 3 c 4 t2 + 4%4 g g3fg2 

7C 1 

°/c4 ft? 2 c 3 ft? 7°/c4c 3 tl 2 ft? 47c4c 3 tlft? 2 16%4tlc 3 ft? 2 
°/cl + °/cl °/cl + 7cl 

_ i3%4^ _ ^um3 _ 62c4 w3 tf _ 26c4 tf + 1Q t3 . 2 ci t22 

/cl /cl 
5°/c4c 3 ft? 3 4_%4cft? 5 c /c4cft? c / c 4 ctl 2%4c 3 ftf 3 4 3 

°/cl + %1 " °/cl + °/cl + °/cl -^ U » cW 
+ 70 c 4 tl 2 ft? 2 -12 c 4 t2 3 tl + 16 c 4 tl 2 t'2 2 + 5 c 4 ft? 4 - 18 c 2 ft? 2 + 2 c 2 tl 2 
+ 46c 4 tl 2 ft? ft? + 4c 4 tl t2 2 ft? - 50c 4 tl t2 ft? 2 + 10c 2 ft? 2 + 22c 4 tl 4 + 5c 4 ft? 4 )/( 

4 %1)) + c 2 tl 2 - 2 c 2 tl ft? + c 2 t2 2 )W) + ^- + \{tl- ft?) c(4min(-(50 

+ 26 c 2 tl ft? - 8 ft? c 4 ft? 3 - 30 tl c 2 ft? + 10 ft? c 2 ft? + 10 ft? 3 c 4 ft? 

6%3c 3 tl ft? ft? 4 7c3 ft? c 3 ft? 2 %3 ft? 2 c 3 ftg 7%3c 3 ft? 2 ftg 

7cl + %1 %1 + °/cl 

4 7c3 c 3 tl t'2 2 16 7c3 ft! c 3 ft? 2 13 °/c3 c 3 tl 2 ft? m , 1 , „ c /c3 c tl 
ft- 1 fr- t^- 62c 4 tl 3 ft? + 



7cl 



2 7c3 c 3 tl 3 5 %3 c 3 ft? 3 4 3 4 %3 eft? 5 %3 c ft? 



7cl %1 %1 %1 

+ 10 ft? 2 c 4 ft? 2 - 30 tl c 4 ft? 3 + 70 c 4 tl 2 ft? 2 - 12 c 4 ft? 3 tl + 16 c 4 tl 2 ft? 2 + 5 c 4 ft? 4 

- 18 c 2 ft? 2 + 2 c 2 tl 2 + 46 c 4 tl 2 ft? ft? + 4 c 4 tl ft? 2 ft? - 50 c 4 tl t'2 t3 2 + 10 c 2 ft? 2 
+ 22 c 4 tl 4 + 5 c 4 ft? 4 ) /(4 7cl), -(50 + 26 c 2 tl ft? - 8 t3 c 4 ft? 3 - 30 tl c 2 ft? 

+ 10 ft? c 2 « + 10 ft? 3 c 4 ft? + 4%4 gf^ - %U f/ t2 + 7%4 g" 2 ^ 

/cl Tel /cl 

4%4c 3 tl ft? 2 16 %4 tl c 3 ft? 2 13 %4 c 3 tl 2 ft? 6 %4 c 3 tl ft? ft? 
7cl + 7cl %l %l 

- 62c 4 tl 3 ft? - 26c 4 tl 3 t'2 + 10 ft? 2 c 4 t2 2 - 5% ^^ + ^ - ^ 

7cl Tel /cl 

+ + 2%4 f tl3 - 30 tl c 4 ft? 3 + 70 c 4 tl 2 ft? 2 - 12 c 4 ft? 3 tl + 16 c 4 tl 2 t'2 2 

/cl /cl 

+ 5 c 4 ft? 4 - 18 c 2 t'2 2 + 2 c 2 tl 2 + 46 c 4 tl 2 ft? ft? + 4 c 4 tl ft? 2 t,? - 50 c 4 tl ft? ft? 2 

+ 10c 2 ft? 2 + 22c 4 tl 4 + 5c 4 ft? 4 )/(4 7cl)) + c 2 tl 2 -2c 2 tlt2 + c 2 ft? 2 / 1 / 2 ) 

7cl :=2c 2 tl 2 + oc? ft? 2 + 5 c 2 ft? 2 - 9 - 8 ft? c 2 ft? - 2 tl c 2 ft? - 2 c 2 tl ft? 

7c2 := 50 + 20 c 2 tl ft? - 10 ft? c 4 ft? 3 - 20 tl c 2 ft? + 50 ft? c 2 ft? + 10 ft? 3 c 4 ft? - 56 c 4 tl 3 ft? 

- 16 c 4 tl 3 t,? + 8 ft? 2 c 4 ft? 2 - 30 tl c 4 ft? 3 + 68 c 4 tl 2 tS 2 - 10 c 4 ft? 3 tl + 8 c 4 tl 2 ft? 2 
+ 5 c 4 ft? 4 - 35 c 2 t,? 2 + 32 c 4 tl 2 ft? fcg + 14 c 4 tl t£ 2 ft? - 46 c 4 tl t'2 t3 2 -15 c 2 ft? 2 
+ 18 c 4 tl 4 + 5 c 4 ft? 4 - 4 c° tl 0 - 34 c° tl 4 ft? ft? + 34 c° tl 3 ft? ft? 2 - 14 c° tl 2 t3 3 t2 

- 30 c e tl 2 ft? 2 ft? 2 +34 c° tl 3 ft? 2 ft? - 14 c° tl 2 ft? 3 ft? + 6 c° tl 3 ft? 3 + 12 c c tl 5 ft? 

- c e tl 2 ft? 4 - 13 c e tl 4 ft? 2 +6 c° tl 3 ft? 3 + 12 c° tl 5 ft? - c° tl 2 ft? 4 - ft? 4 c G t2 2 

- 2 ft? 3 c e ft? 3 - ft? 2 c e ft? 4 - 13c 6 tl 4 ft? 2 + 10c 6 tl t? 3 ft? 2 + 2c c tl ft? 4 ft? 
+ 10 c G tl ft? 2 ft? 3 + 2 tl c e ft? 4 ft? 

7c3 := 26 c 3 tl 2 ft? - 14 c 3 tl 2 t2 + S c 3 tl ft? 2 - 32 tl c 3 ft? 2 + 2 ft? 2 c 3 ft? -8 1,? c 3 t^ 2 

- 4c 3 tl 3 + 10c 3 ft? 3 + 10c ft? -8ct2- 2ctl + 12c 3 tl t^ ft? + 6 V%2 

7c4 := 26 c 3 tl 2 ft? - 14 c 3 tl 2 t2 + 8 c 3 tl t^ 2 - 32 tl c 3 t3 2 + 2 ft? 2 c 3 t^ - 8 ft? c 3 t^ 2 
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41. Graphical Methods For Target Position Report And Error Likelihood Ellipse Display 

The analytical methods demonstrated here relate the graphical depiction of the Error Likelihood Ellipse to 
the relative position of a given Target Position Report and the most likely position value within the Error 
Likelihood Ellipse. These depictions are not meant to list an exhaustive summary of all such graphical 
methods. The policy of the PQIC OSMMTS System is to leave such methods to the particular needs of an 
implementation's circumstances. The methods found herein are for demonstration purposes only. 

> with(plots) : 

Warning, the name changecoords has been redefined 

> F : =PL0T3D (POLYGONS ( seq ( seq ( [ 

> [(i-l/2)/mesh, (j -1/2) /mesh, MZ[1] [i,j]] , 

> [(i+l/2)/mesh, (j -1/2) /mesh, MZ[1] [i,j]] , 

> [(i+l/2)/mesh, (j+l/2)/mesh,MZ[l] [i,j]] , 

> [(i-l/2)/mesh, (j+l/2)/mesh,MZ[l] [i,j]]] , 

> i=MZ[2] . .MZ[3]),j=MZ[4] . .MZ[5])), 

> COLOR (ZHUE) , AXESSTYLE (NORMAL) , AXESLABELS("X" , "Y" , "Z") , 

> VIE¥(TPR[l]-l/2. .TPR[l]+l/2,TPR[2]-l/2. .TPR [2] +1/2, DEFAULT)) : 

> SI : =PL0T3D (POINTS ( [MZ [6 . . 7] /mesh , 0] ) ,C0L0R(RGB ,0.0,0.0,1.0), SYMBOL (CR0 

> SS,50)): 

> S2 : =PL0T3D (POINTS ( [TPR [1 . . 2] , 0] ) , COLOR (RGB ,1.0,0.0,0.0), SYMBOL (CROSS , 5 

> 0)): 

> display ({F, SI, S2}); 




> F:=PL0T3D (POLYGONS ( seq ( seq ([ 

> [(i-l/2)/mesh, (j-l/2)/mesh,MZ[l] [i-1, j-1]] , 

> [(i+l/2)/mesh, (j-l/2)/mesh,MZ[l] [i,j-l]] , 

> [(i+l/2)/mesh, (j+l/2)/mesh,MZ[l] [i, j]] , 

> [(i-l/2)/mesh, (j+l/2)/mesh,MZ[l] [i-l,j]]] , 

> i=MZ[2]+l. .MZ[3]) ,j=MZ[4]+l. .MZ[5])), 

> COLOR(ZHUE) , AXESSTYLE (NORMAL) , AXESLABELS ( "X" , "Y" , "Z") , 

> VIE¥(TPR[l]-l/2. . TPR [1] +1/2, TPR [2] -1/2. .TPR [2] +1/2, DEFAULT)) : 

> S : =PL0T3D (POINTS ( [TPR [1] , TPR [2] , TPR [3] ] ) , COLOR (RGB ,1.0,0.0,0.0), SYMBOL 

> (BOX, 50)): 

> display ({F,S}); 
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(which would be an "infinite-many sided polygon"), and the n-many SDU's would be at positions 

2k^i ( (2k-K\ . . /2A-7T 

re « — r cos +tsin 

V V n J \ n 

for k — 0, 1, 2, .... n — 1, and where r is determined by the effective range of the possible target positions, 
relative to an arbitrarily chosen "origin" for location purposes only. Note how such an arrangement au- 
tomatically satisfies Placement Rules #1 and #3, and Placement Rule #5 is violated only when a target 
comes between the SDU's on the "edges" of the containing circle. This consideration may be addressed 
through the choice for r, which also covers Placement Rule #2. 

The value of r is called the Radius of Coverage, and the placement of the units so described is called 
the SDU Arrangement 

For example, for n — 8, the SDU's would be at positions e Z x , for k — 0, 1,2, . . . , 7. This translates to 

r (cos 77 + i sin ?/) 

where 

7T 7T 37T 57T 37T 7?T 

"- U T2- l"- 1'2 - 1 

(see Figure 13). 

Note how, for n — 8, there is a problem with having SDU's at both 



^cos ^ + i sin ^ and r 



37T . . 3tt 
cos — — + 1 sin — 
4 4 



since sin | — sin ^ , which violates Placement Rule #4. However, choosing any prime n > 4 satisfies 10 
Placement Rule #4. 

The optimal placement of SDU's may now be characterized analytically in terms of the "point of 
diminishing returns." 

Condition 24 Given a coverage function A(n.r) that describes the area covered by the n-many SDU ar- 
rangement with radius of coverage r. find 

(OA 

it* (;• I s c ) - mm < — (n. r) < s c 



Then given r, (r) is the minimum number of SDU's necessary to satisfy the placement rules, and 
where the coverage area instantaneous rate of increase first drops below £ units of area. By choosing a 
sufficiently robust value for £, the optimal number 71* may be explicitly calculated. 



56. Formulation of A (n, r) 

Since the SDU arrangement is completely symmetric for any two SDU's and the center of the containing 
circle, consider the sector spanned by two SDU's A and D, and the center of the containing circle O (see 
Figure 14). Note that the position and orientation of this figure is arbitrary, and the results of this section 
do not depend on the particulars of such aspects. 

Lemma 25 The area of AAOB is \r 2 sin (w (l - £)). 



and 



& 2 K 
da 2 



^ < (vs+ « )<0 

a=b=c 16 (&0?f U ' 1 3 ^ V5/ 



regardless of the value of a. Similar calculations confirm a! % < 0 and , -i\ < 0 regardless of the value of b or 

00 \a=b=c oc \a=b=c 

c. 

15 If n is prime, then 

2u>l7T , 2w-2TV 
Wl 7^ W-2 => COS ^ COS 

71 71 

and 

, . 2u.'l7T , . 2u.'27T 

wi 7^ W2 =>- sm 7£ sin 

n n 

since n would not be divisible by 2. or by any w-i. Hence, no two SDIPs would have the same x-valuc or same y-valuc. 
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2. The OSMMTS System Interface 

The OSMMTS System Interface consists of one Principal Application Specific Integrated Circuit (PASIC) 
Central Processing Unit (CPU), generically referred to as the PCPU, a set (at least four) of remote sensing 
SDU's that send information to the PCPU, and a database of statically stored data that the PCPU accesses 
for parameter 3 data, algorithm exceptions, and other information, which are used to produce the Target 
Position Report (TPR) 4 , as well as supporting reports as the implementation determines (see Figure 2). The 
PCPU, SDU's, and any database systems must be coordinated on and agree with an absolutely maintained 
time system, accurate to at least twice the precision of the anticipated Target Position Report. In this 
respect, the accuracy of the Target Position Report depends more on the maintenance of the time system 
than it does on the algorithmic introduction of error due to calculation round-off or lack of precision. 

A Target Position Report is generated whenever a SDU sends a coordinated stream 0 of timing infor- 
mation to the PCPU. Since different SDU will send information at slightly different times about the same 
target, an absolute timing schedule must be used to ensure valid comparison of timing data from the SDU 
set. 

A Target T may only initiate a signal to the SDU set when t — 0 mod where <f — 10 n / p cycles in a 
10 n Hz PCPU, where there are p signals per second For example, for a 1 GHz PCPU, if a target sends a 
signal to the SDU set every h, second, then p — 2, and £ — — w l^ L0 2 — 10 9_1 ° Sl0 2 . 

The Effective Range of the OSMMTS System is the maximum time for this receive/query/confirm 
period. It measures the farthest a target may be away from the closest qualifying set of SDU's and still be 
detected by the system. 

A complete Signal Period, i.e., £ — 10 n /p cycles in a 10™ Hz PCPU, consists of six Phases, each 
encompassing an interaction between the PCPU, the SDU set, and the parameter database (see Figure 1). 

1. Receive, during which the PCPU receives the detected signal information from the SDU set. This phase 
must last as long as the effective range, plus overhead time for communications between the SDU set 
and the PCPU. The information passed during this phase consists of: 

(a) SDU ID 

(b) Target ID 

(c) Time Of Signal Detection 

The SDU and Target ID are static codes used throughout all phases and signal periods. If either 
the SDU ID or the Target ID changes during a signal period, it must be through a formal change 
management process incorporated into the particular implementation of the OSMMTS System. It shall 
be the responsibility of the OSMMTS System implementation to ensure that changing SDU ID and/or 
Target ID are linked properly for inference purposes. 

The Time of Signal Detection is relative to the common absolute timing mechanisms in the OSMMTS 
System. 

2. Query, during which the PCPU queries the sending SDU for a confirmation code to ensure communica- 
tion integrity. If the confirmation code sent by the SDU is not correct (see the next phase), the PCPU 
queries the SDU again for the proper confirmation code. This is repeated up to a tunable number of 
iterations. If no correct confirmation code is received in the allotted time, the SDU is deactivated. 

3. Confirm, during which the PCPU receives and processes the confirmation code sent by the SDU. It is 
during this phase that any required re- transmissions are also requested, received, and disposed. 

4. Process, during which all calculations are completed to produce the Target Position Report, and sub- 
sequent reports for evaluation, quantification, and adjustment purposes. 

0. Report, during which the Target Position Report and supporting information are made available on 
output channels, and during which any auxiliary communications with the SDU's is completed. 

6. Sync, during which no processing activity is scheduled. This is useful when coordinated processing 
activities require synchronized signal periods. 

One signal period begins when the previous one ends. The sync phase may be used to coordinate 
any overhead processing issues to implement this requirement. 

3 Parameter data consists of Demerit, History, and Confirmation values for each SDU. Sec Section 45 on page 38 for more 
information. 

4 Sec Section 42 on page 37. 

* In the OSMMTS, "coordinated" means "timed to a cycle coordinated by a master clock." 
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6. Calculation Algorithm 

To summarize, given SDU positions (*i,yi), (x2-,y2)-, and (x 3 ,y 3 ), and observed data {t1.t2.t3}, the following 
algorithms give the Calculated Target Position (xo,yo) of Target T. 



7. Common Elements 

A. Calculate fx — t\ — t 2 and f 2 — tx — t 3 . 

B. Calculate dx-J (x L - x 2 ) 2 + (yx - y 2 ) 2 and d 2 - J (x-x - x 3 ) 2 + (yx - y 3 ) 2 . 



8. When y x + y 2 , Vi + 2/3, and d t > \ fi\ c, for i - 1,2 

1. Calculate mi = -(^). 

2. Calculate m 2 = - (f^-). 

tse Oi - 2(3,2-3,!) 

. T . , /2c((/ 1 +/ 2 )c+V , 4 K +(/ 1 c) 2 )+( 3; l- 3; f)+(3,|-3,;) 
4 ' tSe ° 2 - 2(3,3-3,1) 

5. Solve 

= «ff - (/l c) 2 / (-1 - (^)) (*> - (^Efc)) 

2 1 + (* - (-1 (^fe) + M) (* - h {^k) + *)) 

for /t; call this value k 0 . 

6. Evaluate 61 and b 2 with k — kq: call these values /?i and /? 2 , respectively. 

7. Then 

(*o-.2/o) = ; mi 



9. When 1/1 - 1/2, 2/i ^ 2/3, and > |/;| c, for i - 1,2 



1. Use 7"i — — — 



/icV4« 



2(x 2 -xi) 

/ 



2. Use r 2 - 

3. Solve 



2c ((fx + / 2 ) c ± v^7(Mfj + (4 - x\) + (yg - yf) 
- (±fxC^AK + (fx^ 2 + (4-4)) (S) 

rf! - (Ac) 2 



w ~~ ~ 2 + (^l - 7 "i) («2 - ri) + (1/1 - r 2 ) (y 2 - r 2 ) 
for k; call this value Ko. 

4. Evaluate ri and r 2 with k. — kq-. call these values 71 and 72, respectively. 

5. Then 

(xq-.Vo) - (71.72) 



26 
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-4c 3 ti 3 + 10c 3 t3 3 + Wct3 - Sct2 - 2ctl + 12c 3 tl t2 13 -dV%2 



27. Assembly Calculation Sequence For A Target Position Report 

The order of arithmetic steps to calculate a Target Position Report, given arrival times ti,t2,t;3, does not 
depend on the coefficients found in a particular MAPLE calculation. In fact, a generic calculation sequence is 
available regardless of the values assigned to the (xi, r/i). For the purposes of this section only, the coefficient 
and pi values for the X -value of Target Position Report Example #1 shall be used as a demonstration of the 
calculation sequence that may be utilized to assemble a Target Position Report. This sequence is oriented 
to the anticipated MMIX assembly language implementation that forms the PQIC OSMMTS instantiation. 



°/cl := oc 2 t2 2 + 2c 2 tl 2 + oc 2 t3 2 - 9 - 2 tl <? 13 - 2c 2 tl t2 - 8 t3<? t2 

%2 := 50 - c 6 tl 2 t2 4 + 12 c 6 tl 5 t2 + 6 c° tl 3 t2 3 - 13c 6 tl 4 12 2 + 12 c G tl 5 t3 - c 6 tl 2 t3 4 
+ 6c c tl 3 t3 3 - 13 c G tl 4 t3 2 - 16c 4 tl 3 t2 - 14c 6 tl 2 t3 3 t2 + 34c 6 tl 3 t3 t2 2 

- 34c 6 tl 4 13 t2 - t3 2 c G t2 4 - 2 13 3 c° t2 3 - t3 4 c G t2 2 + 34c 6 tl 3 t3 2 t2 

- 30 tl c 4 t3 3 + 8 c 4 tl 2 t2 2 - 30 c e ti 2 t3 2 t2 2 - 20 ti c 2 - 10 c 4 t2 3 tl 

- 14 c° tl 2 t2 3 t3 + 20 c 2 ti ti? + 10 c G ti t2 2 t,? 3 + 2 c 6 ti t2 4 13 + 10 c G t2 3 t3 2 
+ 68 c 4 tl 2 t3 2 + 2tl c G t3 4 ti? + 8 t3 2 c 4 ti? 2 +d0t3c 2 t2- 10 *3 c 4 ti? 3 - 4 c G ti G 

- 56 c 4 ti 3 t3 + 10 t3 3 c 4 ti? + 5 c 4 t3 4 - 35 c 2 1^ 2 + 32 c 4 ti 2 t5 t# + 14 c 4 ti ti? 2 t3 

- 46 c 4 ti t2 t3 2 -15 c 2 t,? 2 + 18 c 4 ti 4 + 5 c 4 t2 4 

%3 := 26 c 3 ti 2 tS - 14 c 3 tl 2 t2 + 8 c 3 ti t^ 2 - 32 ti c 3 t3 2 + 2 t^ 2 c 3 t2 - 8 i3 c 3 ti? 2 

- 8cti? + 10c 3 t,? 3 - 2cti + 10 ct3 + 12c 3 ti t3 -Ac 3 tl 3 - 6 V%2 



%4 : = %3 + 12vYc2 




+ 16 c 4 ti 2 ti? 2 - 30 ti c 2 t,? - 12 c 4 12 3 tl + 26 c 2 ti ft? + 70 c 4 ti 2 t,? 2 




.2 



4%4c 3 ti ti? 2 7%4c 3 ti 2 tj? „ 0 , ^ 3 xo 6%4c 3 ti fc? 
%1 + %1 %1 



+ 10 i3 3 c 4 fc? 
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r0.025 

^0.02 

r0.015 z 

r0.01 
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B 
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o 



Figure 14: Optimal Placement Sector 



Proof. In Figure 15, ACOB is a right triangle with height OC and base CB. Since OB — OA. then 
/? — 7r ~ Q , corresponding to Z.OAC and Z.OBC. 
Then 

0(7 — rsin./? and CB — rcos/? 

which means 

Area Of AC05 = 2 r 2 sin/? cos/? 
Since ACOB is similar to ACCM, then 

Area Of AAOB = r 2 sin./? cos./? = 2 ?- 2 sin 2/? 

In the SDU arrangement, we have a — 2 £ . for n > 4. Hence, the area covered by the convex hull of 



SDU's is incrementally given by * ?* 2 sin 2./?, where 



7T- ( 2 f) _ n7r-27r _ / n-2 \_ A _ 1 
2 2n ~ V 2n ) ~ * [ 2 n 



Hence, for an n-many SDU arrangement, we have 

Area Of AAOB = ^r 2 sin 28 = ^r 2 sin (n (l - 



Claim 26 A (n, r) = ' 2 V 2 sin (tt (l - 2 )) 

Proof. The coverage region is comprised on n-many exactly sized sectors of the type found in Figure 
14. The result follows immediately by Lemma 25. 
For a fixed r, as n oc, we have 

2^ 



TT 1- 



and 



sin 7r 1 — 
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Signal Timing 




Figure 1: The OSMMTS Signal Timing Cycle 



A special SDU, the Calibrating SDU (CSDU), is placed where the PCPU can always receive its signal 
as clearly as possible, i.e., without corruption. This CSDU calibrates the system by being a known distance 
from all other SDl/s, and allows for a consistent, guaranteed error-free TPR by which all other SDl7s are 
evaluated for accuracy and consistency. 



16 
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10. When y x ^ y 2 , yi - y$, and (k > |/^| c, for i - 1, 2 

1 T'«» - /2C 2 (/i+/2)±V4*o+(/ic) 2 +(z!-^) 



, , ±/i<V4«o + (/ic) 2 + (ai-af) + (^-^) 

2- Ue r 2 - 2{y2 _ yi) y_^,.^ f ../,,, . ;;; ; ) 

3. Solve 

« - -J — + (»i - 'l) (»2 - 'l) + (yi - r 2 ) (2/2 - r 2 ) 
for re; call this value kq. 

4. Evaluate ri and r 2 with re — reo; call these values 71 and 72, respectively. 
0. Then 

(#0,1/0) = (Ti.Ti) 



11. Calculated Target Position Examples 

Suppose (#1,1/1) — (2,3), (#2,1/2) — (4,2), and (#3,1/3) — (3,1), and observed data t± — 1.6, £2 — 2, 
£3 — 1.75. For demonstration purposes only, let c — 1 (see Figure ). Then 1/1 7^ 1/2, Vi 7^ 1/3, and > c, 
fori = 1,2. 

1. We have /1 = t x - 1 2 = 1.6 - 2 = -0.40, and / 2 = ti - t 3 = 1-6 - 1.75 - -0.15. 

2. Then d x = ^(2 - 4) 2 + (3 - 2) 2 = Vo and d 2 = ^{2- 3) 2 + (3 - l) 2 - yo. 

3. And mi = - (^if.) =2. 

4. Then wis = - (fl 2 ) - \. 

_ .. . -0.40A/4re+(-0.40) 2 + (4 2 -2 2 ) + (2 2 -3 2 ) o 7 , , . . . 

0. Also 61 — 2(2-3) 25 v2ore + l — 2 . Aote that the positive sign has been 

chosen. 

-0.1o((-0.40-0.15)+v'4«;+(-0.40) 2 ) + (3 2 -2 2 )+(l 2 -3 2 ) ? 

6. And fr 2 2(1 _ 3) '— — - + J0V20K+ 1. Aote that the posi- 
tive sign has been chosen. 

7. Solving for re 0 gives us 

8. Therefore 



= 26454174801 561 V 292162046321 = 1.5369135 • 

18790926400 2348865800 



2 /27205811857 561 7 



-a — 



+ _„„V292162046321 - = -2.9976991 • • • 



and 



9. So 



25 V 751637056 93954632 2 

1167 3 /27205811857 561 , nnn ^ nnAc .^ no ,„.-„ 
^ = 1600 + 200 V 751637056 + 93954632 ^292162046321 = 0.823oo64 •• • 



(# 0 . vo) - ( H ~ • my ( h ~ h \ ■ h) - 1(2.5475037 • • • , 2.0973083 • • • )| 
to 7 decimal places (see Figure 8). 

Now suppose (#1,1/1) — (2,3), (#2,1/2) — (4,3), and (#3,1/3) — (3, 1), and the observed data is ti — 1.6, 
t 2 — 2, t 3 — 1.75. Note that 1/1 = y 2 and 1/1 ^ 1/3, and > |/;| c, for i — 1, 2. For demonstration purposes 
only, let c — 1. 

1. We have /1 = t x - 1 2 = 1.6 - 2 = -0.40, and f 2 =h~ t 3 = 1.6 - 1.75 = -0.15. 

2. Then d x = y/{2- 4) 2 + (3 - 3) 2 = 2 and d 2 - \/(2 - 3) 2 + (3 - l) 2 = Vo. 

3. This means ri — _0 - 40 V 4re +(-°-40) +(ie-4) _ _ i_^2ok + 1 + 3. Note that the positive sign has been 
chosen. 
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4 %4 ft? c? t'2 2 . 4 4 _ 5 %4 c ft? %4 c « _ 5 %4 c 3 ft? 3 4 %4 c ft? 
+ %1 + ° C %1 + %1 %1 %1 

2 %4 c 3 tl 3 

+ 0/1 -18c 2 ti? 2 + 2c 2 tl 2 + 46c i tl 2 t3t2 + 4c 4 tlt2 2 t3 

/cl 

- 50 c 4 tl tj? ft? 2 + 10 c 2 ft? 2 + 22 c 4 « 4 + 5 c 4 ft? 4 ) /(%1) 



#2 := -(50-26c 4 tl 3 t£ 

- 30 tl c 4 ft? 3 + 16c 4 tl 2 t2 2 - 30 tl c 2 ft? - 12c 4 ft? 3 tl + 26c 2 tl fc? + 70c 4 tl 2 t3 2 
+ 10 13 2 c 4 ft? 2 + 10 13 c 2 t2 - 8 t3 c 4 t£ 3 - 62 c 4 tl 3 + 10 13 3 c 4 ft? 

_ 4 °/c3 c 3 tl ft? 2 7 %3 c 3 tl 2 fcg _ 13 %3 c 3 tl 2 ft? _ 6 %3 c 3 tl ft? tS 

%1 + c /cl " %1 ~ %1 

16 %3 tl c 3 ft? 2 4 °/c3 ft? c 3 ti? 2 _ %3 ft? 2 c 3 ft? . 4 4 %3 c tl 

+ °/cl %1 %1 + ° C W + %1 

5%3c 3 ft? 3 4%3cft? 2°/c3c 3 tl 3 5%3cft? 1D „ o2 „ 

- %1 + %1 + %1 - %1 -^^^+2c 2 tl 2 

+ 46c 4 tl 2 ft? t2 + 4c 4 tl t2 2 ft? - 50c 4 tl t2 13 2 + 10c 2 ft? 2 + 22c 4 ft/ 4 + 5c 4 i? 4 )/( 
%1) 



#3 := c 2 tl 2 - 2 c 2 tl ft? + c 2 ti? 2 
#0 := (2 tl-tg- ft?) c 



#00 :=-*(« -tf)c 



A := (#0 + V#l + #2 + #3) x #00 



#5 :=|(«-tg)c 



#6 : = - (50 - 13 %A t i " 2 * 5 " 26 c 4 tl 3 tg - 30 tl c 4 ft? 3 + 16 c 4 tl 2 ti> 2 - 30 tl c 2 
/cl 

- 12 c 4 t£ 3 tl + 26 c 2 tl *g + 70 c 4 tl 2 ft? 2 + 10 t3 2 c 4 t'2 2 + 10 ft? c 2 t2 - 8 *3 c 4 ft? 3 
°/ c 4 ft? 2 c 3 ti? 16 %4 tic 3 ft? 2 4 C / C 4c 3 tlt£ 2 7%4c?tl 2 t2 nn 4i< ^ 0 

— %r~ + — %i %r— + 62c4 " 

_ 6 %4 c^tl ft. « + 1Q ^ c4 tf + 4 %4 c 3 t. 2 + . c4 ^ _ 5 %4 c t3 
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42. The Target Position Report Vector 

The Target Position Report Vector (TPR) is the quadruple defined by the Calculated Target Position 
(xo,yo), the Maximum Likelihood Target Position (x' 0 ,y' 0 ), associated error likelihood ellipse + 

(m°I™ 2 ) 2 = and the associated likelihood value for {mi ^^ R ~ M 3 (ft,1%,t$): 

tprv = , (f:sj + (* :»)* = * <mi ,;r^ - * «•«.«) 

The value of i? is chosen so that (mi ~™ 2)2 -R is at the 100(1- a) % percentile of the M 3 (tl,t? 2 ,t* 3 ) 
distribution. 
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In fact. 



Figure 15: Sector Division 



lim 'V 2 sin (it (l — 2 )) 

n — i. rv~- ^ * ^ ft// 



which is the area of a perfect circle of radius r. 
Now by Claim 26, we have 



n— >-oc 

— r 2 hm 



)) 



7t— ^tx; „j 

- "^Mi™ COs(7r(l- £)) 

- -Trr 2 (-1) 



— 7rr* 



fj(n,r)=r 2 Q sin (vr (l - £ 



H — COS T 1 

n \ V tt. 



the graph of which is contained in Figure 16. This may be used to calculate n* for a given £. 
Finally, solving A (n, r) for r gives 



r — \ I A esc ( 7r ( 1 — 
/ n V V n 



and if the approximate coverage area A and the approximate radius of coverage r are given, then define 



71 „ 



,-Acsc//- 
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SRI In 




Figure 2: The OSMMTS Processing System 



3. Target Position Calculations 

4. Introduction 

Suppose there were two Surface 0 Detection Units (SDU) labeled I and II that detect a signal that propagates 
through a uniformly dense environmental medium 7 encompassing the SDU and a Target T 8 . The signal 
from the target is detected at times t± and t-z- where the detection is made under a common methodology 9 . 
Suppose further that the SDU are a known distance d > 0 apart, and that T and the SDU are .?i > 0 and 
s 2 > 0 distances apart, respectively. Finally, let 0 be the angle (measured positively) between T and the 
SDU (see Figure 5). In the case that T is on the line segment between I and II, take 9 — | , i.e., the triangle 
formed by I, II, and T is degenerate, with two 0 angles between I and T, and II and T. If I, II, and T are 
otherwise collinear, then take 8 — 0: the other angles will be irrelevant. 

The calculations that follow will depend on axioms that collectively assert the OSMMTS System. These 
axioms formally state the conditions or assumptions under which the OSMMTS System may be implemented. 

Axiom 1 Given d, and the positions of I and II, the only observed data, and therefore the only basis for 
any calculations found in the OSMMTS, are the times ti and t- 2 of signal detection. 



6 The detection units arc called "surface" units because all coordinates under the OSMMTS may be thought of as (x, y, z) € 
R 3 , were z — 0. Even after rotations and translations, the coordinates mav be thought of as planes or surfaces embedded in 
R 3 . 

7 The term "environmental medium" is intentionally generic. The analytical and implementation aspects of the OSMMTS 
apply equally well to signals broadcast through the air, under water, in gas-filled containers, in vacuums, or in any other 
medium. 

However, the medium will be assumed to be of uniform density to ensure the validity of the formula: 

Distance — (Linear Speed) * Time 

The "uniformly dense" restriction may be relaxed in subsequent OSMMTS System versions. 

8 Even though this development applies to a single Target T, all analytical methods in the OSMMTS may be extended to 
any number of individual targets with distinguishable signals. 

9 The "common methodology" refers to the electrical, mechanical, physical, or other well-defined, reproducible, and deter- 
ministic means of declaring a detection complete at an SDU. Whatever particular methods arc used, the implementation must 
be uniform and constantly applied to each SDU in the OSMMTS System to ensure the integrity of algorithmic calculation. 
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Figure 8: Example Target Calculation Result when yi ^ y-z- yi ^ y?,-, and d{ > \fi\ c. for \ 



4. And 



r-2 



' -0.15 ^(-0.40 - 0.15) ± y/4n + (-0.40) 2 ^) + (9 - 4) + (1 - 9) 

- (-0.4Oy/4n + (-0.40) 2 + (16 - 4)) ( 
| (-0.15 * 0.55 - O.Oov^k + 0.16 + 9) 

i mn — inn V^O/C + 1 



1 

2(1-3) 



— 1C00 200 

Note that the positive sign has been chosen. 
5. Solving for k 0 gives us 

5958301793 1233^23857734 

K 0 - 



6. Therefore 



and 



3769960000 117811250 
1 / 6109100193 1233^23857734 



71 25 V 150798400 H 



3567 



72 1600 200 V 150798400 H 



4712450 
6109100193 1233^23857734 



+ 3 = 2.7414202 • 



4712450 



= 2.1975504 • 



7. So 



(a; 0; yo ) = (2.7414202 • • • , 2.1975504 • • • ) 
to 7 decimal places (see Figure 9). 



12. Algorithm Implementation Summary 

( TPR Calculation 16 ) 

( TPR Example #1 SDU Coordinates 18 ) 

( TPR Example #2 SDU Coordinates 23 ) 

( ELE Calculation 34 ) 

(ELE Example #1 SDU Coordinates 36) 

( ELE Example #2 SDU Coordinates 39 ) 

( Mitigations For Reflections 50 ) 

( Coverage Area Optimization 58 ) 
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%4c« o%4c 3 t3 3 4%Act2 2%4c 3 tl 3 1Q „ o2 n _ - 

+ ^%1 %1~ + ^1~ + ^1 ^« 2 +2C« 2 

+ 46c 4 tl 2 t3 t'2 + 4c 4 tl t» 2 t3 - 50c 4 tl t£ tS 2 + 10c 2 13 2 + 22c 4 tl 4 + 5c 4 t£ 4 )/( 
%1) 



#7 := -(50 - 26 c 4 tl 3 ftg - 30 tl c 4 tS 3 + 16 c 4 tl 2 t'2 2 - 30 tl c 2 tS -12 c 4 t? 3 tl 
+ 26 <? tl t2 + 70 c 4 2 t? 2 + 10 tS 2 c 4 t? 2 + 10 tS c 2 t£ - 8 tS c 4 ft? 3 - 62 c 4 tl 3 

4 ^ 4%3c 3 tlt? 2 1%Z<?tl 2 t2 \Z%Zc?tl 2 t3 
+ ™*** t *- %l + %i " % i 

_ 6 %3 c 3 tl ft? t3 16 %3 tl c 3 t3 2 4 %3 *3 c 3 t'2 2 _ %3 t,? 2 c 3 t2 . 4 4 
°/cl + °/cl %1 %1 + ° C W 

%3 c tl _ 5 7c3 c 3 tS 3 4 °/c3 c fcg 2 %3 c 3 tl 3 _ 5 %3 c t3 2 2 

+ %1 %1 °/cl °/cl %1 

+ 2 c 2 tl 2 + 46 c 4 tl 2 t5 t2 + 4 c 4 tl t2 2 t5 - 50 c 4 tl t£ t,? 2 + 10 c 2 t? 2 + 22 c 4 tl 4 

+ oc 4 fc? 4 )/(%l) 



#8 := c 2 tl 2 - 2 c 2 tl «0 + c 2 ti? 2 
B:= (#5xV#6 + #7 + #8)+#4 
x 0 :— A + D 
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43. Containment Policies 

44. Definitions and Remarks 

Definition 14 A TPR is said to be accurate */ the actual position of the target is inside the associated 
ELE. The value of R determines the likelihood of this event. 

Definition 15 Any calculation algorithm used to produce a set of numerical values intermediate and inferior 
to the TPR is called an analytical step. 

Definition 16 An analytical step is called a mitigation if it is taken before the time data {h,t-2, . . . ,tk, . . .} 
are collected. 

Definition 17 An analytical step is called a optimization */ it occurs after the time data {h,t-2, ...,tk,.. .} 
are collected. 

Remark 18 The purpose of mitigation steps is to reduce the error variance a. 

Remark 19 The purpose of optimization steps is to increase the likelihood of an accurate TPR. 

Definition 20 An irregularly occurring, non-analytical step taken at any time to accomplish the same goals 
as mitigation and optimization is called ad-hoc. 

Definition 21 The collection of ad-hoc, mitigation, or optimization steps taken in an implementation of the 
OSMMTS is called the system 's containment policies, and referred to individually as a system containment 
policy. 

45. The Demerit System 

The OSMMTS Demerit System is an ad-hoc containment policy that acts simultaneously as a mitigation 
and an optimization. Under this system, the three SDU's chosen to calculate the TPR are those three that 
are most likely to produce the "best" TPR based on past performance (thereby making it an optimization 
step), by way of reducing the variability of the utihzed data (thereby making it a mitigation step). 

Suppose there are n-many SDU's, however, only k < n many receive a signal within the reception 
window. There are (^)-many combinations of SDU's, and (^)-many combinations of the A;- many that 
receive the signal taken three at a time. Each SDU has three values associated with it at the beginning of 
each processing cycle, namely its non-negative Demerit Count, its positive History Total, and its possibly 
null Boolean Confirmation Value. At the beginning of all processing, the demerit count for each SDU will 
be zero, the history total with be one, and the confirmation value will be null. The confirmation value at 
the beginning of the processing cycle is determined by its observed value during the confirmation cycle. At 
the end of a processing cycle, the demerit count and history total are determined by the steps below, and 
the confirmation value is set back to null. 

For each processing cycle, and for each of the (g)-many combinations, the following steps determine the 
end-of-processing-cycle demerit counts and history totals. 

1. Set the likelihood value A. 

2. Eliminate those ri-many combinations that are colhnear. 

3. Eliminate those T2-many combinations that do not all have positive history totals and TRUE confirma- 
tion values. 

The SDU's involved in the (ti + r2)-many combinations eliminated in Steps 2-3 are called deficient 
for the current processing cycle. This designation is removed at the beginning of a new processing cycle. 

4. Among the remaining, i.e.. qualifying combinations, choose the combination of three SDU that collec- 
tively have the minimal sum of demerits. 

5. In case of a tie in Step 4, use the combination with the largest history sum. In case of a further tie, 
choose the combination with the smallest individual demerit count. In case of a last tie, randomly 
choose uniformly among the finalists. 

The combination so chosen is called the calculating combination, and the SDU's involved are called 
the elected SDU's. Increment the history total by 1 for each elected SDU. 

6. Subtract two demerits from the count for each elected SDU. Recall the demerit count for an SDU cannot 
become negative. 
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5. Derivation 

A critical element of the OSMMTS is to calculate a Target Position Report, which contains a position report 
based on the arrival times at the SDU's. Consider arrival times ti and t- 2 at two SDU, named I and II, for 
the same Target T. 




Figure 3: Relative placement of I, II, and T 



For the triangle formed by I, II, and T depicted in Figure 5, the Law of Cosines gives 

d 2 - sf + s% - 2sis 2 cos 9 (1) 

Let / — t\ — 1-2, c — the speed with which the signal linearly propagates through the environmental 
medium, and let t be the time at which the Target T sends the signal, with units in common. Note how 
ti—t only when T is identically collocated with either I or II. The term "speed" refers to the instantaneous 
change in distance per unit time, and it is expressed in terms of two-dimensional Euclidean geometry. 

Axiom 2 For the purposes of the OSMMTS. it will be assumed, without loss of generality, that ti > t. for 
t = l,2. 



Then 

or 

Hence, 



fc - {h-t 2 )c 

- ((h-t)-(t 2 -t))c 

- (h -t)c- (t 2 ~t)c 

- §1 - S 2 

s'f - 2s lS2 + sj - (fc) 2 

s'f + s'l - 2sis 2 cos# - d 2 — 0 
s{ + s 2 - 2*!* 2 - (fcf - 0 
2s-i 5- 2 (1 - cos6>) ^ d 2 - (fc) 2 

d 2 - (fcf 



(2) 
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Figure 9: Example Target Calculation Result when 1/1 — y-i. y± ^ j/3, and di > \fi\ c. for % — \. 2. 



13. Target Position Report Calculation 

A Target Position Report is the calculated (x, y) position value based on the arrival times detected at the 
SDU's, and the position information for the SDU's involved in the calculation. There are no adjustments 
made for the uncertainty in the arrival time data, as such issues are addressed by the Error Likelihood 
Ellipse methods. 

The principal calculation policy in the PQIC OSMMTS System is to base a Target Position Report on 
the choice of SDU arrival times that maximizes the likelihood that the calculated Target Position Report 
is as accurate, i.e., close to the actual target position, as possible. Since the PQIC OSMMTS SDU's do 
not change position as a function of time when the system is operating, this optimization may be achieved 
by implementing the PQIC OSMMTS Target Position Report Standard Methodology, which means pre- 
establishing the calculation algorithm that the system will utilize when a particular combination of SDU 
arrival times are used to calculate the Target Position Report. 

When error conditions are calculated, e.g., in the Error Likelihood Ellipse, the same algorithm set shall 
be used with appropriately adjusted arrival time data. 



14. Target Position Report Standard Methodology 

The PQIC OSMMTS Target Position Report Standard Methodology for calculating a Target Position Report 
begins with the use of a position-specific combination of the SDU's taken three at a time. Based on this 
position information, the Target Position Report algorithm calculates a (x, y) position value for the given 
*i»*2»*3 values. This algorithm is specific to the combination of SDU's that receive the arrival time data. 
Therefore, for n-many SDU positions, there must be (3) -many algorithms pre-established for use by the 
Target Position Report calculation device. For example, if there are 5 SDU positions, then ( 3 J 10 ) separate 
algorithms must be established to calculate a Target Position Report regardless of which three SDU's are 
involved in the calculation. 

The set of algorithms generated for a particular set of SDU positions differ only by the number of terms 
and the coefficients of the terms involved in the sum. An effective organizing system that addresses the 
numerical issues related to such a complicated sum may be found later in this document. 
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28. Error Distribution 

All data collected during process actions are subject to error. The source or cause of the error varies by 
implementation, however, the OSMMTS, like all such processes, has several dominant sources of data error. 
There is mechanical error introduced by the physical limitations of the mechanical devices to precisely 
calculate the exact time of signal detection. There is propagation error from the natural corruption of the 
signal before detection. There is also rounding error due to the finite precision calculations utilized in the 
implementation, and algorithm error introduced by inefficient numerical methods. 

The PQIC OSMMTS may be used to specifically addresses rounding and algorithm errors by framing the 
implementation of all algorithms in assembly code, e.g., within a system that uses the MMIX 11 instruction set 
and memory model, even though this approach is not required in any particular OSMMTS implementation 12 . 
However, utilizing MMIX-based assembly code may be considered the OSMMTS "preferred embodiment," 
with regard to algorithm implementation, since all non-mechanical and non-propagation error issues are 
addressed by this approach. 

Finally, propagation error may be addressed, if not eliminated, by the optimal placement and mitigations 
for reflections methods documented in this memorandum. Mechanical errors may be addressed by the 
merchantability of the PCPU. 

The primary focus of the error distribution calculations is to account for the remaining effect of error 
in the arrival time data, regardless of source, after all precautions, optimizations, and mitigations have been 
used to minimize said error. 

Definition 8 The MLTarget Position 'is the calculated target position most likely to be correct given the 
arrival time data and the variability of the arrival times. It is found at the same time the M 0 distribution 
is calculated (see 9). 

29. Derivation 

Suppose ti and t-2 are not precisely measured, in the sense that 

t{ —ti+ ci and *2 — t-2 + £-2 
where the t| are the observed values, and the U are the actual "true" times. Here, each Si is an error term. 
Suppose each Si ~ N (0, a 2 ) . for a common a > 0, and all such e; are IID. Then 

/* - i* - 1% 

= (h+e 1 )-(t 2 +£ 2 ) 
= {h-t 2 ) + {si-s 2 ) 

= f + (E 1 -E 2 ) 

which means 

/* -/ = c!-c 2 ~iV(0,2(7 2 ) 

Then 

s* — s* 2 — f*c and s L — s 2 — fc 

mean 

(«i - 4) - (*i - *2) - (/* - /) o ~ N (0, 2cV 2 ) 

or 

(*2i " *ii) " («2i - *n) - - (A* - /i) c - N (0, 2cV) 

(*22 - *12) - (*22 - «2l) - - (f*2 ~ H) C ~ N (0, 2cV) 

Now from (1), we have 

df — s 2 ! + s 21 - 25n §21 cos(9 0 

and 

df = «Ii+«2i -2siiS^cos^ 

11 MMIX (2009): A "RISC Computer For Tlic Third Millennium, as developed by Dr. Donald E.~Kmith of Stanford University, 
is documented at http://www-cs-facTil.ty. Stanford. edu/'kmith/mmix. html. 

12 An OSMMTS implementation will work as well without MMIX as any other implemented OSMMTS as long as the error 
issues identified herein arc addressed. 
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since 9 ^- 2nir, for any integer n. If I, II, and T are collinear, e.g., the Target T is on the line between I 
and II, then take 9 — it (as in the case of a degenerate triangle). Assume otherwise 10 that I, II, and T are 
not collinear. 
Note how 

2(1- cos 9) 2k 

Here, 

re > 0, when d > \ f\ c 
k — 0. when d — |/| c 
re < 0, when d < |/| c 
since d and c are both necessarily positive. 
Now $i — $2 — fc gives us 

s{ ~ SCSI ~ « 
2 



Si - 

Hence. 




si - ^ ^/c ±\Jak + (/c) 2 ^ and s 2 - 



T > (4) 

The choice of sign in (4) depends on the signs of / and re: whichever makes si > 0 and s 2 > 0 should 
be used. 

However, since si > 0 and s 2 > 0, then re ^ 0, hence 

2 K _ ( 2k \ ( JcTV'4 K +(/c) 2 

f 2k \ ( fcTV^+ifc) 2 

2/t(/cTV , 4/t+(/c) 2 ) 
— =4^ 

(/c T VW(/c) 2 ) 
2 



So 



- i(-/c±/l«+(/c) 2 ) 

Sl - 2 (/<= ± \/ 4 * + (/ c ) 2 ) 311(1 6 '2 = 2 ( _/c ± V 4k + C/ c ) 2 ) ( 5 ) 

Expression (5) shows that the are functions of re, /, and c. Since re is a function of d, /, c, and #, 
we have the as functions of rf, /, c, and 9. The values of d and c are fixed by the circumstances and by 
nature, and the / value is observed; this means the Si may be thought of as functions of 9 given d and the 
data ti and t 2 . 



We have 

4« + (/<) 2 = 4 (1^) + (/c) 



2 

' rf 2 -(/c) 2 \ , (l-cos<?)(/c) 2 
, l-cos6> y ' 1-cosO 
-2(/c) 2 + (l-co^)(/c) 2 
1-cos 9 
2</ 2 -(l+cosg)(/c) 2 
1— cos# 



hence 



10 If I, II, and T arc collinear, yet the Target T is not on the line between I and II, then while it would make sense to take 
0 — 0, this would make n ill-defined. For the purposes of the OSMMTS development, assume I, II, and T arc collinear only 
when the Target T is on the line between I and II. 
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15. Target Position Report MAPLE Algorithm 

Given the position information of the SDL"s that receive the arrival time data, the following MAPLE algo- 
rithm calculates the Target Position Report based on the arrival times. For the purposes of this document, 
only example algorithms shall be given based on example position information. The completely generic 
form of the Target Position Report MAPLE Algorithm shall be documented in a separate PQIC OSMMTS 
technical memorandum. 

In general, each term of the Target Position Report calculation has the form 

(sign) x (coefficient) x (ratio) x <f° x t L pl x t 2 P2 x h Ps 
for pi — 0, 1, 2, .... where po — pi + p2 + P3- The sign is either positive or negative 1, the coefficient is a 
positive integer, and the ratio is either 1 or a ration of previously calculated subexpressions. If the coefficient 
is allowed to be nonnegative, i.e., take on the value 0, then, in general, each Target Position Report may be 
expressed as the sum of the same number of such terms, where the number of such terms depends on the 
SDU position information. 

The Target Position Report terms may be organized into convenient subassemblies for calculation 
efficiency by using the tables found below as worksheets, as demonstrated later in this document. 

16. { TPR Calculation 16 ) = 
osmmts : — proc (nl ,n2, n3 ) 

local fl ,f2 ,dl ,ml , m2 ,bl,b2, k, kOx. betal . beta2 : 

global c, xO ,y0 .xl ,yl ,x2 ,y2 , x3 , y3 ,sl ,s2,s3,k0, cthO , thO : 

options ' Copyright 2003PQI Consulting All Rights Reserved': 

(/e5cnpfrwn"0SMMTSuTargetuPositionuReport u anduSupportinguInforiiiation": 

fl: — nl — n2: f2: — nl — n3: 

dl : - sqrt((x2 - xl) © 2 + (y2 - yl ) © 2): 

ml: = -((x2 - xl)/(y2 - yl)): 

m2 : = -((x3 - xl)/{y3 - yl)); 

bl: ={jl * c * sqrt(4 *k + {fl * c) © 2) +{x2 © 2 - xl © 2) + {y2 © 2 - yl © 2))/(2 * {y2 - yl )); 

b2: - (J2 *c*((fl +f2)*c+sqrt(4*k+{fl *c)©2)) +{x3 ®2-xl ®2) + {y3 ®2-yl ©2))/(2*(f/3 -yl)): 

Wx: = ° 

{solve(k= {(dl ©2 -{fl *c)©2)/2) +(xl - ((b2 - bl) / (ml - m2))) *(x2 - ((b2 - bl)/(ml -m2))) 
+(yl - (ml *((b2 - bl)/(ml -m2)) + bl)) *(y2 - (ml *((b2 - bl)/(ml - m2)) + bl)),k) 0 }: 

if (nops(k0x) — 0) then 

k0: — undefined: 

elif (nops(k0x) > 1) then 

if (Mr [2] — conjugate(kOx[l])) then 

k0: — undefined: 

else 

k0\ — min(op(l,k0x),op(2,k0x)); 

end if : 

else 

kO: = k0x[l]: 
end if : 

betal: - eval(bl,k = kO): 
beta2: — eval(b2 ,k — kO); 
xO: — (beta2 - betal) /(ml - m2); 
y0: — ml * xO + betal : 

cthO: = 1 - ((dl © 2 — (fl * c) © 2)/(2 * kO)): 
thO: — arccos (cthO): 

si : - sqrt((xl -x0)®2 + (yl - y0) © 2); 
s2: = sqrt((x2 - x0)(&2 + (y2 -y0)®2): 
s3: ^ sqrt((x3 -x0)®2 + (y3 - yO) © 2): 
return xO ,y0 ,k0 , kOx , nops (kOx). cthO . thO .si .s2.s3: 
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Sl {9\d,{h,h})-\ Uh-t 2 )e± } /^^^-^ >0 

s 2 (o\d, { h , t 2 }) = i ht 2 -t l) c± v ^-^r-S"^ ) > o 



Note how 



Hence. 



s i s 2 — 4 



* (/ciy^+Cfc) 2 ) -^-/ c ±^/4« + (/c) 2 y 

! / (/c) 2 ± 2/c^4« + (/c) 2 + 4k + (/c) 2 
4 \ - (-/c) 2 ± 2fcy/4K + (/c) 2 - 4k - (/c) 2 
- ±/c V / 4k+(/c) 2 



2 2 

1 z Si— S2 



±/cvW(/c) 2 

/ /c 

- ±^4k+(/c) 2 



(6) 



\ ^ ^ ^ 
\ 

si\ ^mll (x 2 ,V2) 




Figure 4: Intersection of Circles d and C 2 

Since si and «2 are distances from T to I and II, respectively, then T must be at the intersection of 
two circles centered at I and II (see Figure 4). In coordinate geometry, suppose I is at {xi,y{), II is at 
(*2»W2)i and T is at (x, y). The origin for the coordinate system is arbitrary, and its position may be chosen 
for convenience in any particular implementation. Then 

(x - x±f + {y- Vl f = sf 
(x - x 2 f + (y - y 2 f = 4 



must be satisfied simultaneously for a given 9. Hence, from (6), 
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8000 
5000 




Figure 16: in-Many Vertex Convex Area Of Radius r 

to be the optimal number of SDU's, where z — fi satisfies 16 

A (z — 7r) + irr 2 sin z — 0 



57. MAPLE Implementation Code 

> X : = (n/2) *r ~2*sin (Pi* ( 1- (2/n) ) ) ; 



> solve (A=X , r) [1] ; 



> subs(A=1000,n=9,A=X) ; 

> f solve C/.,r,0. .100); 



1 . 2 
X :— nr 2 sin(7r (1 — )) 

2 n 



n 



1000- !J r 2 sin( 77r ) 



18.59345062 



16 Note that if A (fi — ir) + wr' 2 sin //, — 0., then 



2 _ -§Acsc/ 



?r 2 sin(7r(l- ^))| 2 . 

2 \ \ 77. / / I n.— A esc /i 



7'" sin I 1 



r A esc /T. 



= Acsc/z.sin^l--^)) 

— A esc //. sin ( ( 7r — -jj- 7rr 2 sin //,) ) 

= A esc // sin ((7r + (//, — 7r))) 

= A esc // sin fi. 

= A 
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7. Calculate the TPR using the calculating combination. 

8. Calculate the A-ELE for the calculating combination. 

9. Calculate the ^(*) — (ti + r 2 )^-many TPR for all other qualifying combinations. Each of these TPR 
is called an Alternate Position Report (APR). 

10. For each APR calculated in Step 7, if the APR falls outside the A-ELE, then add one demerit to the 
count for each SDU involved in the APR. 

11. For each APR calculated in Step 7, if the APR falls inside or on the A-ELE, then subtract one demerit 
to the count for each SDU involved in the APR. Recall the demerit count for an SDU cannot become 
negative. 

12. Add one demerit for each SDU that does not report a positive confirmation. 

13. When the demerit count for an SDU exceeds the Warning Threshold, send an alert to report a frequently 
deficient SDU. 

14. When the demerit count for an SDU exceeds the Terminal Threshold, shut down communication with 
the SDU and do not consider it further (by setting its history total to zero) until explicitly reset. Also 
send an alert to report a failed SDU. 

15. These steps are in addition to the disabling of an SDU if proper query responses are not confirmed 
during the receive phase. 
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> solve (A=X,n) ; 

2A 

r 2 sin(RootOf(A .Z + irr 2 sin(_Z) - Ait)) 

> subs(A=1000,r=18.25,A=X); 

1000 = 166.5312500 nsin(7r(l - -)) 

n 

> fsolve(7„,n,4. .100); 

12.10617840 

> plot3d(X,r=0. ,10,n=3. .20, 

> orientation=[55,60] ,axes=normal,grid=[10,20] 

> ); 



r300 
^250 
-200 
M50 




> dXdn:=diff (X,n); 

2 

1 2 r 2 cos(7r(l ))ir 

dXdn := - r 2 sinfyr (1 - -)) + ^ 

2 n n 

> plot3d(dXdn,r=0. .100,n=3. .50, 

> orientation=[25,60] ,axes=noraal,grid=[10 

> ,20]); 



> C:=subs(r=25,dXdn) ; 

r , 625 . , „ 2„ 625co 5 (^(l-^))^ 

C := — sm(7r (1 - - + 

2 n n 

> solve (C=10, n) ; 

-2 7r/(RootOf((15625sin(_2') 2 - 15625) _Z 2 + (-31250 tt sin(_Z) 2 + 31250 tt) _Z 
+ 15625 tt 2 sin(_^) 2 - 15625 tt 2 + 15625 sin(_^) 2 - 1000sin(-Z) + 16) - tt) 
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so 



However, 



s ll + s 21 



and 



Therefore, 



so 



5 11 + S 21 ~ 2 «o cos 9q — «n + «21 _ 2k o cos ^0 

|(/ic±\^+(i^)) + (± (-/iC±^4«o + (/ic) : 

= i(4(/ lC ) 2 + 8«o) 
- (/ic) 2 + 2« 0 

*u + 4 2 = (/i*<) 2 + 2kS 

(Ac) 2 + 2« 0 (1 - cos0 o ) - {flcf + 2k* (1 - cos^) 

1 



k* (1 - cosflj) - * 0 (1 - cos^o) - 2 (((A* - & - s ^ c ) 2 - (A*c) 2 ) (12) 
These calculations lead to the first important subsidiary distribution for quantifying the OSMMTS 
error. 

Definition 9 Let rj — rj 01,02,03) — cos#q — cos# 0 - Then v is said to have an Mo (t*,t 2^3) 

distribution. 



Hence, 
and 



1 — cos# 0 — 1 — cos0q + n 



(*2? - - (4i - 4i) = T/rc^S + (/fc) 2 ± Acy/iKo + (Ac) 2 

= ta*V 4 *o + arc) 2 ± ac^2 (;tio J :£) + m 2 



— < 



±(A*-(o!-02))c 



2 / ^-((/r-( gl - g2 )) C )- 

\^ (l-COS^S+T/) 

M +((/r-(ci-c 2 ))cr 



If & = (A* - (01 - 02)) c ~ iV (A*c, 2c 2 o- 2 ), then we have 

0 M-(A*c) 2N 



(^-*I?)-(*ii-*n)==FA*c 



(A*cr±?iW2 



^(l-COB^) 

Definition 10 WTien £ ~ A< (A*c, 2c 2 ct 2 ) andi] ~ M 0 (tf,q,t^). then 



(l-cos6£+r ? ) 



(13) 



2 (2/2 - 2/i ) 



cos ^ + rf) 



■a 



is distributed according to the Mi (t*,t2>*si) probability distribution function. 

We also have that when £ 2 = (A* ~ ( £ i ~ ^?)) c ~ N (/2* c > 2cV 2 ) and rj ~ M 0 {tl,t* 2 ,tl), then 



Definition 11 If X ~ My {f,\.t* 2 .t*^) andY ~ Mi (**,£■!, tlj), and X andY are independent, then 

x-y~ M 2 (ti,r 2 ;t* 3 ) 

Definition 12 // J ~ M 2 {tl.t* 2 .t%) and K ~ M 2 (t*,t 2 '*3)> ar ^ ^ and K are independent, then 

( mi _ m2 ) J2 + ^1 _ ^ 2 ^ K 2 ~ Af3 
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46. Mitigations For Reflections 

Given a TPR (xo,yo) and multiple detection times at an SDU, say {t1.t2-.t3; ■ ■ .}, where t\ is used to generate 
(xo, 1/0)5 the question of mitigating the error in the TPR due to reflections 1 ' 3 of the "true" signal may be 
stated as follows: What is the earliest time a reflection of the "true" signal could be detected at an SDU, 
and what is the latest such time. Any signal detected between these two extremes would be considered a 
reflection, and therefore not used in the calculation of a TPR. Any signal detected "too early" or "too late" 
would indicate the detection time data is corrupt and should not be used to calculate the TPR. 

Axiom 22 // the signal detection times recorded during the Receive Phase are inconsistent with a calculated 
TPR, the TPR should not be used. 



47. Definitions 




Figure 10: Generic Straight Line Obstruction 



Consider a generic straight line obstruction as depicted in Figure 10. Even though the following develop- 
ment requires a straight line obstruction, a more generalized parameterized formulation may be substituted 
for 14, such as a Bezier curve or trigonometric function. 

As s goes from 0 to 1, point (x s , y s ) goes from L — (xi, yi) to R — (x r , y r ), so that 

(x s ,y s ) - s(x r ,y r ) + (l-s)(xi,yi) ^ 
- {xi + s{x r -x{),yi + s{y r -y{)) K ' 

Then given SDU position (xi,y L ). TPR position (xo,yo). and obstruction endpoints L and R, we wish 
to find 

dmm - min {di + d 2 } 

0<s<l 

and 

<W - max {di + d 2 } 

0<s<l 

subject to 

d 2 -d\+dl- 2dxd 2 cos (ir - 2a) 

13 The term "reflections" also includes signal behavior characterized by "cchos," "secondary," and "phantom." 
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58. (Coverage Area Optimization 58} = /* Embedded In MAPLE Code*/ 
This code is used in chunk 12. 



59. Sample Calculations And Results 

> f solve (C=10,n=3. .50) ; 

13.62494565 

> f solve (C=5,n=3. .50) ; 

17.21207895 

> f solve (C=l,n=3. .50) ; 

29.51887628 



